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INTRODUCTION 


ETONATION in the internal combustion engine 
has been studied by many methods since the 
pioneer work of Midgely and Boyd! described in 1922. 
Carbon resistances,”* piezoelectric crystals,‘ condensers,* 
electromagnetic generators,“"* photoelectric cells,’ and 
moving mirrors have been used as sensitive elements to 
pick up the deflections of a diaphragm exposed to 
cylinder pressures. Several types of optical and electri- 
cal systems have successfully produced records from the 
pressure pick-up units. By means of transparent win- 
dows placed in the combustion chamber walls flame 
motion has been recorded,'®''!* combustion tempera- 
tures'? measured, and radiant energy analyzed.'*! 
Flame movement has also been followed by means of ioni- 
zation gaps exposed inside the combustion chamber.'* 

From the detonation data collected by various means 
a consistent theory of the process has gradually evolved. 
This theory holds that detonation, or knock as it is 
often called, is due to a very rapid pressure rise in some 
part of the burning gases which is followed by intense 
pressure waves within the charge.*!”'*'* Most con- 
vincing evidence of this detonation theory is given by 
Withrow and Rassweiler® who discuss results from 
simultaneous flame photographs and pressure records. 
Despite the vast amount of time expended on the 
problem, very little quantitative information has 
appeared on the intensity of the pressure waves which 
accompany detonation and no systematic effort has 
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been made to correlate observed results with charge 
conditions within various combustion chamber shapes. 
The author has devoted part of his time to this subject 
since 1927. A discussion of certain theoretical aspects 
of the problem and experimental results from two en- 
gines were published as an N.A.C.A. Report in 1934." 
This study was handicapped by the lack of an amplifier 
and recording system with uniform sensitivity over a 
sufficiently wide frequency range, while the theoretical 
discussion was not given in a form suited to engineering 
use. Recent work on detonation at the Massachusetts 
Institute of Technology has been directed toward the 
refinement of indicating equipment and the practical 
application of theoretical results. The investigation is 
being continued with grants of funds and equipment 
from the N.A.C.A., the Bureau of Aeronautics of the 
United States Navy, the Sperry Gyroscope Co., and the 
Wright Aeronautical Corporation. 


GENERAL THEORY 


Prediction of the types of pressure waves to be ex- 
pected in the combustion space of an engine after knock 
has occurred requires fitting a solution of the wave 
equation for sound to the given geometrical bounda- 
ries.2!_ This procedure involves assumptions that 
the pressure wave amplitude is small compared to the 
average pressure and that the gas particle velocity 
normal to the walls is zero. The particular case of a 


circular cylinder with flat ends at right angles to the 
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cylinder axis has been worked out by the author and 
applied to the C.F.R. Engine and the N.A.C.A. Uni- 
versal Test Engine. The excellent checks between 
calculated and measured detonation frequencies cer- 
tainly justify further applications of the theory. In 
the present paper the former results will be extended to 
include energy computations and the outline of a 
method for investigating the nature of the initial dis- 
turbance responsible for the pressure waves. 

Based on the earlier work, the instantaneous pressure 
disturbance from equilibrium can be expressed as 


mg=0,1,2... 
mr =0,1,2... 
mg =0.1,2... 
Where 
n, = number of azimuthal nodes (nodal planes con- 


taining the cylinder axis and a diameter). 
number of radial nodes (nodal cylinders). 


n, = 
mn, = number of axial nodes (nodal planes perpen- 
dicular to the cylinder axis). 
Pngn,n, = absolute amplitude of pressure associated with 
the mode n,n,n, (mode with nodal surfaces* 
NN,N,). 
Gngn, = Number associated with the nodal surfaces 
myn, (the method for determining these 
numbers is outlined in reference 19). 
r = radial distance from cylinder axis. 
a = radius of cylinder. 
6 = azimuth angle. 
nn, = phase angle associated with azimuth. 
z = distance from bottom of cylinder measured 
parallel to the cylinder axis. 
h = height of cylinder. 
Mnonn,= frequency of the mode 
= time. 
Engn,n, = phase angle associated with time. 


The frequency 1,, gms is related to the velocity of sound 
c, and a distance \,, wrt, called the wavelength by the 
equation 

= ¢/ Angn,n, (2) 
Sound velocity c can be found from measurements of 
equilibrium pressure Po, equilibrium density pp) and y the 
ratio of specific heat at constant pressure to specific 
heat at constant volume as" 


= V yPo/p (3) 


Values of wavelength in terms of cylinder radius for 
a number of the simpler modes of vibration are given in 
Figs. 1, 2, and 3. 

* As used in the text, the term “nodal surface” refers to a sur- 
face over which no pressure variation due to a particular mode 


occurs. This differs from the usage of reference 19 in which a 
nodal surface refers to a surface of zero particle velocity. 
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Fic. 2. Ratio of wavelength to cylinder radius for 
nm, = 1. 


ENERGY CALCULATIONS 


In general, the energy associated with a system of 
pressure waves in an enclosed volume V is 


W = (1/yP) SSS P?(xy2)dV (4) 


Where p(xyz) is the pressure amplitude considered as 
a function of position only and the integration is ex- 
tended over the entire volume. Applied to the par- 
ticular case of a circular cylinder with flat ends, Eq. (4) 
becomes for each mode of vibration 
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Fic. 3. Ratio of wavelength to cylinder radius for 
nN, = 


It is convenient to express the relationship between 
energy and pressures in terms of a ratio M where 


ngn,n, = Engn,n,/ (pngn,n,/ yP. 0) (6) 


In this equation Ex, n,n, is the average energy density 
associated with the mode 1,n,n,. 

The ratio M is given for a number of modes in 
column 4 of Table 1. Column 5 of this table is the cal- 
culated frequency corresponding to these same modes of 
vibration based on typical conditions existing near top 
center in a C.F.R. Engine with a compression ratio of 
4.85. 


PRESSURE WAVES DUE TO GIVEN INITIAL KNocK 


The type and amplitude of the vibrations produced 
by a given knock will depend upon the shape, location, 
and intensity of the initial disturbance. For computa- 
tion it is assumed that this disturbance can be approxi- 
mated by assigning a constant pressure p) to a region 
between 7, and a, 0 and 2, —6, and @,. If both sides 
of Eq. (1) are multiplied by the variables on the right- 
hand side and integration over the whole cylinder 
volume is carried out for ¢ = 0, an expression for the 
pressure amplitude of each mode of vibration is ob- 
tained as 


1 21a 


n/a 


TABLE 1 
Based on Sound Velocity of 3000 ft./sec. 
in C.F.R. Engine. 
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(5) 


Azimuthal Radial Axial 


Nodes Nodes Nodes Cycles/sec. 
Ng n, n, Magn: Nngnynz 
0 0 0 1.000 13540 
0 0 1 .500 19650 
0 0 2 .500 30830 
0 1 0 . 1620 24900 
0 1 1 .0810 29200 
0 1 2 .0810 37000 
0 2 0 .0903 35900 
0 2 1 .0452 39000 
0 2 2 .0452 46300 
1 0 0 .119 6500 
1 0 1 .0595 15420 
1 0 2 .0595 28900 
1 1 0 .0537 18800 
1 1 1 .0269 23650 
1 1 2 .0269 34200 
1 2 0 .0368 30100 
1 2 1 .0814 34200 
1 2 2 .0814 41900 
2 0 0 .0332 10620 
2 0 1 .0166 17780 
2 0 2 .0166 29650 
2 1 0 .0450 24000 
2 1 1 .0225 28200 
2 1 2 .0225 36400 
2 2 0 .0318 35800 
2 2 1 .0159 39000 
2 2 2 .0159 46300 
3 0 0 14700 
3 1 0 28400 
4 0 0 19000 
4 1 0 33100 
5 0 0 22900 
5 1 0 37000 
6 0 0 26300 
6 1 0 43300 
7 0 0 30000 
& 0 0 33800 
i) 0 0 37700 
bon, [1 a) — 
(7 


Pngnrnz = 


1 1 
0 0 0 
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The integral in the denominator has already been evaluated in connection with the energy factors of Eq (6). 


Using this result the pressure amplitudes can be written 


Po 1 
Prgnyng = J J Ing [angn,(1/a) — (8) 


0 


where V = volume of charge. 


The integrals of Eq. (8) are somewhat difficult to cal- 
culate but the work has already been carried out for 
several cases. Preliminary checks with experimental 
results have shown qualitative agreement but a de- 
tailed discussion of this subject is beyond the scope of 
the present paper. 


EXPERIMENTAL METHOD 


The ultimate object of the M.I.T. Detonation Re- 
search is to produce an instrument suitable for quantita- 
tive measurements on any type of engine under any 
operating conditions as well as in the laboratory. The 
first step in this program was to construct laboratory 
apparatus to record accurately the physical phenomena 
available to operate an instrument in routine work. 
For the present case a reasonably certain check be- 
tween theory and experiment required at least two 
simultaneous records of detonation waves in order to 
measure phase and amplitude relations. With the 
amplitudes of pressure waves small compared to total 
cylinder pressures, it is difficult to study the wave mo- 
tion: directly from a pressure card. Because of the 
emphasis on the pressure waves found in a rate of change 
of pressure record,**:!® the author chose to continue the 
development of an electromagnetic instrument which 
gives an output directly proportional to the rate of 
change of pressure. 

In order to insure that no essential part of the detona- 
tion process is distorted the indicator must give true 
records of pressure variations from those of the ordinary 
indicator card up to frequencies of about 25,000 cycles 
per sec. This requires a diaphragm natural frequency 
in the vicinity of 100,000 cycles per sec. in order to pre- 
vent errors from mechanical effects.* In considering 
errors of this type the diaphragm can be considered as 
substantially undamped. A second advantage of using 
a high natural frequency for the diaphragm is that the 
effects of engine vibration are minimized. 


Pick-Up UNITs 


Fig. 4 shows the essential parts of two types of elec- 
tromagnetic pick-up units. In both cases the necessary 
magnetic field is supplied by a central core of heat- 
treated cobalt steel. The right-hand unit has an annu- 
lar air gap with a cylindrical inductor machined as an 
integral part of the diaphragm. Any deflection will 
cause a transfer of flux from one side to the other of the 
stationary coil and thus produce an induced voltage. 
The left-hand unit also has a stationary generating coil 


Fic. 4. Electromagnetic indicators for rate of change 
of pressure. 


wound on the central core and the flux changes are 
due to changes in the air gap between the pole piece 
and the diaphragm in the manner of the conventional 
telephone receiver. A water jacket around the magnet 
is shown in both units. Many hours of operation in 
various engines have shown that the water jacket is un- 
necessary in practice. The permanent magnet seems 
to retain its strength satisfactorily after an initial 
“shakedown’’ period. 

It is possible to construct a very sensitive unit using 
a cylindrical inductor but the natural frequency is of 
necessity much lower than that possible with the flat 
diaphragm type. This natural frequency advantage 
makes the telephone receiver type more desirable than 
the other instrument. The diaphragm in use at the 
present time is 0.060 in. in thickness and 0.500 in. in 
diameter which gives a natural frequency of about 
95,000 cycles per sec. Photographs of the actual in- 
strument are shown in Fig. 5. 


AMPLIFIER-OSCILLOGRAPH SYSTEM 


The wide frequency range required makes it necessary 
to use cathode ray oscillographs for the indicating ele- 
ments. For satisfactory photographic records with the 
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CIRCUIT DIAGRAM OF HIGH GAIN AMPLIFIER 


Fic. 6. Wiring diagram of high gain section of amplifier. 


are 
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CIRCUIT DIAGRAM OF HIGH VOLTAGE POWER SUPPLY 


Fic. 7. Wiring diagram of high voltage amplifier and 
power supply. 


film speed of 400 in. per sec. a specially designed power 
supply of 4000 volts was used. The amplifiers were 
specially designed in two sections, a high-gain section 
connected between the indicators and a high-voltage 
section to operate the cathode ray tube. Figs. 6 and 7 
show wiring diagrams of the amplifier units. Sensitivity 
curves for the two units are given in Fig. 8. 


WHEN RELATIVE AMPLITUDE =! 
INPUT VOLTAGE |2.5 MILLIVOLTS 
ATTENUAITION 
AMALIFIER! | 
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© +3 Au | 
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APPLIED FREQUENCY CYCLES PER SEC 
Fic. 8. Overall sensitivity curves of amplifier—oscillo- 
graph systems. 


CAMERA 


CATHODE Ray SL) 
OSCILLOGRAPH 


CFR. ENGINE 


Fic. 9. Arrangement of apparatus for study of detona- 
tion in the C.F.R. Engine. 


Fic. 10. Photograph of apparatus for recording detona- 
tion in the C.F.R. Engine. 


ENGINE TEST 


The general arrangement of apparatus for tests on a 
C.F.R. Engine is indicated in Fig. 9. A special spacer 
between the cylinder and the head was used to permit 
the location of indicators and spark plugs in four al- 
ternate positions. Fig. 10 shows the actual appearance 
of the setup. 
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Fic. 5. Photograph of rate of change of pressure indicator. : 
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Indicator calibration by non-firing operation in 


Fie. 11. 
the C.F.R. Engine. 


CALIBRATION 


Overall calibrations of the apparatus are carried out 
by comparing indicator records from the apparatus with 
arate of change card derived from an M.I.T. Pressure 
Indicator record.* Fig. 11 shows the result of such a 
calibration. By this method the maximum sensitivity 
was found to be 4.3 X 10-4 mm. per lb. per sq.in. per 
sec. This value can be varied over a range of one 
thousand to one by means of the attenuators. 


TYPICAL RECORDS 


Fig. 12 shows records from successive cycles in the 
C.F.R. Engine with moderate detonation as judged by 
ear. Tests under identical operating conditions with 
blind plugs shielding the indicators from cylinder gases 
showed that direct vibration effects were negligible. 
Records taken at the same engine speed but without 
audible detonation were substantially free from vibra- 
tion during the power stroke. 


OBSERVED FREQUENCIES 


The most prominent frequency shown in Fig. 12 is 
near 6500 cycles per sec. and decreases as the piston 
descends. This value agrees with results already pub- 
lished by the writer and other investigators.5*"9 By 
reference to the last column of Table 1, which is based 
upon a sound velocity corresponding approximately to 
top center power stroke conditions in the C.F.R. 
Engine, the mode of vibration excited most strongly 
has one azimuthal node and no axial or radial nodes. 
Fig. 13 shows the lines of constant pressure amplitude 
for waves of this type. Dotted lines on one side of the 
nodal plane indicate the same time phase while full 
lines on the other side of the node show the opposite 
phase. It is apparent that indicators located at posi- 


tions A, and A, will show substantially equal amplitudes 
of opposite phase so far as this particular mode of 
vibration is concerned. Indicators located at B, and B, 
will show a strong vibration and no vibration, respec- 
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Fic. 12. Typical records of two successive cycles in the 


C.F.R. Engine. Attenuation factor = 200. Speed 

1230 r.p.m., spark advance 39°; inlet temperature 150°F.; 

72 octane fuel; moderate detonation induced by ethyl 

nitrite; mixture ratio 14:1; air consumption 1.13 X 1073 
Ibs. per cycle. 


tively, from this same mode. Thus the upper record 
of Fig. 12 which shows very little of the lowest frequency 
mode at unit 2 can be explained on the assumption that 
one end of the pressure node was located immediately 
in front of this indicator. The lower record of Fig. 12 
in which the vibrations are out of phase corresponds to 
the case with indicators located symmetrically on either 
side of the nodal plane. 

High-frequency vibrations are present on all the 
records especially in the initial stages. An extended 
examination of many records from various engines has 
shown one vibration component of about 50,000 cycles 
per sec. frequency which can possibly be attributed to 
longitudinal vibrations of the indicator pole piece. 
However, there are definitely other frequencies present 
corresponding to more complicated modes of vibration 
as shown in Table 1. Quantitative work on these modes 
would require an increase in film speed of the recording 
apparatus. 


PRESSURE WAVE ENERGY 


One of the most serious difficulties in measuring the 
physical quantities associated with detonation is due 
to the relatively great variations between successive 
cycles even under carefully adjusted engine conditions. 
Until engine operation can be more precisely controlled 
any estimate of the energy associated with detonation 
pressure waves must be considered as only somewhat 
better than an order of magnitude rather than a precise 
value. With this in mind the wave energy associated 
with the cycle recorded in the lower traces of Fig. 12 
will be calculated. 

To determine the pressure amplitude associated with 
each component taken from the records, it is necessary 
to divide the measured rate of change amplitude by 27 
times the corresponding frequency. This means that 
the actual pressures associated with high-frequency 
components will be small unless the record amplitudes 
are large. In the present case only the pressure waves 
of the lowest frequency in Fig. 12 will make any im- 
portant contribution to the energy. Estimating from 
the record that the maximum rate of change amplitude 
near 16° from top center is 2 mm., for the lowest fre- 
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Fic. 13. Lines of equal pressure amplitude for the lowest 
frequency mode of vibration in a circular cylinder. 


quency mode, the pressure amplitude will be about 23 
Ibs. per sq.in. The measured pressure for this crank 
angle was about 350 Ibs. per sq.in. and the temperature 
about 4000°F. absolute as taken from the Hottel chart.*4 
This corresponds to a ratio of specific heats of about 
1.24. From Table 1 the value of Myo is 0.119. 
Substituting these values in Eq. (6) and solving for 
Ej gives an average energy density associated with the 
pressure waves of 1.7 ft. Ibs. per cu.ft. In the C.F.R. 
Engine at 10° after top center this corresponds to a 
total energy of 0.01 ft. Ibs. (1.3 X 10-5 B.t.u.) The 
measured weight of air in the theoretical mixture per 
cycle was 1.13 X 10~* lbs. which gives a total of about 
1.38 B.t.u. chemical energy per cycle*‘ or an available 
energy of 0.35 B.t.u. if a thermal efficiency of 25 per- 
cent is assumed. This indicates that less than 0.01 
percent of the available charge energy is involved in 
the pressure waves accompanying detonation. This 
does not mean that only this small fraction of the 
charge explodes suddenly in the initial disturbance 
since the efficiency of converting chemical energy into 
sound energy is probably low. E. S. Taylor has given 
an estimate for a theoretical cycle with instantaneous 
combustion of the detonating part of the charge limiting 
the highest possible efficiency for sound wave production 
at 13 percent, for a fuel-air ratio of .078 and a compres- 
sion ratio of 6. No doubt combustion is sufficiently 
slow so that the true value is much less. Even after a 
factor of this magnitude is applied, the amount of 
charge required to produce observed detonation effects 
will be a fraction of one percent. 


DISCUSSION AND CONCLUSIONS 


The general experimental conclusion that the energy 
associated with pressure waves due to detonation is a 
small percentage of the total energy falls nicely into 
agreement with the known result that the output of a 
well-designed engine is not affected by moderate deto- 
nation intensities. 


Since the actual wave pressures are too small to di- 
rectly damage engine parts, other effects must be con- 
sidered. Withrow and Rassweiler have advanced the 
theory that in detonation ““..... The increased jacket 
temperature and loss of power both undoubtedly result from 
increased heat transfer from gases to the walls. It is 
apparent that this heat transfer would be greatly augmented 
by the ‘scrubbing action’ of the high gas velocities induced 
by knock...... 

With the present data it is possible to calculate the 
maximum particle velocity at the cylinder wall due to 
the measured pressure waves. In general the particle 
velocity, v, due to a sinusoidal pressure wave in a gas is*! 


v = —(V p,/po2xn)sin 2rnt (9) 


where p, is the space distribution of pressure amplitude. 

Substituting values for the lowest frequency mode 
as given by the records of Fig. 9 and remembering that 
the velocity normal to the cylinder wall must be zero 
gives a maximum particle velocity of 50 ft. per sec. 
While this figure does not represent a very high speed, 
the oscillatory motion will take place very close to the 
cylinder wall and possibly supply the ‘‘scrubbing ac- 
tion’’ suggested by Withrow and Rassweiler. This 
measured gas velocity due to pressure waves will cer- 
tainly be intensified for surfaces near the openings of 
pockets due to alternating flow in and out with a re- 
sultant heating of spark plug points and similar parts. 
In addition, at the instant of the initial explosion, all the 
energy which is later spread throughout the cylinder 
volume will be concentrated in a relatively small region 
with the result of very high temperatures, densities, 
and surface velocities as pressure waves spread over 
the cylinder. The attendant high rate of heat transfer 
can well be responsible for the local pitting effects ob- 
served in the aluminum combustion chamber surfaces 
of certain engines. 

In regard to the problem of selecting a unit of detona- 
tion intensity suitable for all engines under all condi- 
tions, it seems reasonable to use the average pressure 
wave energy density as the quantitative measure. 
From Eq (5) it is seen that this energy density is pro- 
portional to the square of the pressure wave amplitude. 
Thus if an electromagnetic generator type of indicator 
is used, the associated electrical circuit should first 
integrate to reduce all frequencies to a common basis, 
then “‘rectify’’ and amplify according to a square law 
with a final steady indication proportional to the recti- 
fied output. This general problem is now being studied 
and it is hoped to arrive at a satisfactory answer in the 
not too distant future. 
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Handbook of Aeronautics, Vol. II: Aero Engines, Design and 
Practice, by ANDREW SwaN; Pitman Publishing Corporation, 
New York, 1988; 483 pages, ill., $7.50. 

As in previous editions (1931, 1934) this book consists of a 
more or less complete compilation of descriptive matter, plus a 
considerable amount of statistical data, concerning British aero- 
nautical power plants. Foreign practice is also included, though 
to a much more limited extent. The material is largely obtained 
from published sources, very little original data being presented. 
In addition to chapters on engines and their component parts, 
chapters are included on supercharging, fuel systems, lubrication, 
exhaust systems, heaters, starters, and engine testing. A large 
number of performance curves are presented, both for engines and 
their auxiliaries, together with tables of the ‘‘standard’”’ altitude, 
viscosity conversion charts, physical and thermal properties of 
materials, etc. 

In general, the arrangement is good and the presentation clear, 
which results in a volume of considerable value for reference pur- 
poses. On the other hand, publications of this kind invariably 
suffer from the fact that vital statistics on aviation power plants 
are seldom available until the models in question are virtually 
obsolete. For example, no rating for the Wright ‘‘Cyclone”’ 
in excess of 768 hp. is given, and the table on bearing loads (page 
28) includes such obsolete and unsuccessful examples as the 
U. S. A. W-1, the Packard 2500, etc. Fortunately, the British 
engines used for illustration are somewhat more up-to-date. 

The chief fault of this book, however, is the notable absence 
of critical appraisal of designs and design methods, practice 
both good and bad being presented without comment from the 
author. No advantage is taken of the fundamental principals 
of similitude, and thus very few of the data are reduced to those 
non-dimensional ratios which must be the basis of any intelli- 
gent new design. While this book constitutes a welcome addi- 
tion to the designer’s library, one is reluctantly forced to the con- 
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clusion that a definitive work on aero-engine design has yet to be 


written. 
C. FAYETTE TAYLOR 


Aircraft Propeller Handbook, by Kart HANSSON FALkK; The 
Ronald Press Company, New York, 1937; 99 pages, ill., $4.00. 
The first paragraph of the author’s preface reads as follows: 

“This AIRCRAFT PROPELLER HANDBOOK has been prepared to 
meet the needs of engineers, draughtsmen, students, pilots, ser- 
vice personnel, and others in the aeronautic industry who wish to 
obtain concise, practical data without delving into propeller 
theory.” 

This quotation very adequately describes the Propeller Hand- 
book which Mr. Falk has produced. Rules and formulae for 
determining basic propeller characteristics are given, supple- 
mented by numerous charts, curves, and figures. Starting out 
with the main characteristics necessary from a propeller design 
standpoint, the author follows with concise information on static, 
gyroscopic, and aerodynamic characteristics; a chapter on twist- 
ing moments and strength, and then some very interesting in- 
formation on vibration. He then proceeds with citation of several 
examples to illustrate the uses of the formulae and curves of the 
Handbook, with a final chapter on method of approach for a de- 
signer to use in the selection of a propeller for his airplane. 

This Handbook gives every evidence of being a careful re- 
arrangement of the running notes used by a practical designer, 
which is just what the author is, and therefore it appears certain 
that this small propeller handbook will be found extremely 
useful for the practical man who does not care to delve through 
the more complete textbooks or theoretical treatises on propeller 
design in order to have available concise information for immedi- 


ate use. 
T. P. WRIGHT 


Curtiss-Wright Corporation 
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Boundary Layer in Compressible Fluids 


TH. VON KARMAN and H. S. TSIEN, California Institute of Technology 


Presented at the Aerodynamics Session, Sixth Annual Meeting, I. Ae. S. 
January 26, 1938 


SUMMARY 


The first part of the paper is concerned with the theory of the 
laminar boundary layer in compressible fluids. The known 
solution for incompressible fluids is extended to large Mach’s 
numbers by successive approximation. The compressibility ef- 
fect on surface friction is discussed, and the results applied to esti- 
mate the ratio between wave resistance and frictional drag of 
projectiles and rockets. In the second part the heat transfer 
between a hot fluid and a cool surface, then between a hot body 
and a cool fluid is discussed. A general relation between drag 
and heat transfer as function of Mach’s number is given. The 
limits where cooling becomes illusory because of the heat produced 
by friction are determined. 


HE solution of flow problems in which the density 
is variable is in general very difficult; hence, every 
case in which an exact or ¢ ven an approximate solution 
of the equations of the motion of compressible fluids 
can be obtained has considerable theoretical interest. 
Several authors noticed that the theory of the laminar 
boundary layer can be extended to the case of compres- 
sible fluids moving with arbitrarily high velocities 
without encountering insurmountable mathematical 
difficulties. Busemann! established the equations and 
calculated the velocity profile for one speed ratio. (By 
speed ratio is understood the ratio of the airspeed to 
the velocity of sound.) Frankl? also made an analysis 
of the same problem, however, it is complicated and 
depends on several arbitrary approximations. The sen- 
ior author® obtained a first approximation by a simple 
but apparently not sufficiently exact calculation. 
Hence, in the first part of the present paper, a better 
method for the solution of the problem is developed. 
The boundary layer theory for very high velocities is 
not without practical interest. First, the statement 
can be found often in technical and semi-technical litera- 
ture on rockets and similar high-speed devices that the 
skin friction becomes more and more insignificant at 
high speeds. Of course, it is known that with increasing 
Reynolds Number, the skin friction coefficient is de- 
creasing, 1.e., the skin friction becomes relatively small 
in comparison with the drag produced by wave forma- 
tion or direct shock. Since high-speed flight will be 
performed mostly at high altitude where the air is of 
very low density, so that the kinematic viscosity is 
large, the resulting Reynolds Number will be relatively 
small in spite of the high speed. 
Another interesting point in the theory of the bound- 
ary layer in compressible fluids is the thermodynamic 
aspect of the problem. In the case of lcw speeds the 


influence of the heat produced in the boundary layer 
can be neglected both in the calculation of the drag 
and of the heat transferred to the wall. In the case of 
high speeds, however, the heat produced in the bound- 
ary layer is not negligible, but determines the direction 
of heat flow. In the second part of the paper a few 
simple examples of heat flow through the boundary 
layer are discussed. 

It has been found necessary in most parts of this 
analysis to make the assumption of laminar flow. 
This assumption was found necessary because of the 
lamentable state of knowledge concerning the laws of 
turbulent flow of compressible fluids at high speeds. 
This assumption is somewhat justified by the fact 
that—as mentioned above—in many problems where 
the results of this paper can be applied, the Reynolds 
Number is relatively small, so that a considerable por- 
tion of the boundary layer is probably, de facto, 
laminar. Ackeret‘ called attention to the possibility 
that the stability conditions in supersonic flow might 
be quite different from those occurring in flow with low 
velocities. The authors also believe that the stability 
criteria as developed by Tollmien and others, cannot be 
applied without modification. Finally, some conclu- 
sions of the paper, as will be pointed out, are also ap- 
plicable to turbulent flow. In other cases, as in the 
calculation of drag, the assumption of laminar flow 
surely gives at least the lower limit of its value. 


If the x-axis is taken along the plate in the direction 
of the free stream, the y-axis perpendicular to the plate, 
and « and v indicate the x and y components of the 
velocity at any point, then the simplified equation of 
motion in the boundary layer is 


where both the density p and the viscosity » are vari- 


ables. 
The equation of continuity in this case is 
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A third equation determines the energy balance be- 
tween the heat produced by viscous dissipation and the 
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heat transferred by conduction and convection. With 
the same simplification as used in Eqs. (1) and (2), 
one can write 


_d(,ar ow \? 


where c, is the specific heat at constant pressure, and 
is the coefficient of heat conduction. If Prandtl’s 
number, ¢,u/ is assumed to be equal to 1, then it can 
be easily shown that both Eqs. (1) and (3) can be 
satisfied by equating the temperature 7 to a certain 
parabolic function of the velocity u only. This rela- 
tion between 7 and w is 


To To (= ) U U ( ) 
(4) 

where 

U = free stream velocity. 

M = speed ratio, or Mach’s number of the free 

stream. 
T) = temperature of the free stream. 
T,, = temperature at the wall of the plate. 


Differentiating Eq. (4) one obtains 
To oy w U 2 T, ov w 
(5) 


where the subscript w sefers to conditions existing at the 
surface of the plate. Now (0u/dy),, is always positive; 
therefore, if [(k — 1)/2]M? > (T,,/To) — 1 heat is 
transferred from the fluid to the wall, if [(« — 1)/2]M? 
= (T,,/To) — 1 there is no heat transfer between the 
fluid and the wall, and if [(« — 1)/2]M?< (T,,/To) — 1 
heat is transferred from the wall to the fluid. If there 
is no heat transfer, the energy content per unit mass 
(u?/2) + c,T is constant in the whole region of the 
boundary layer.*® 

The pressure being constant the relation between p 
and T is, 


p= (6) 
The expression for the viscosity based on the kinetic 
theory of gases is 
= (7) 


However, the following formula is in closer agree- 


ment with experimental data 
(7a) 


Busemann! calculated the limiting case for which 
[(x — 1)/2]M? = (T,,/To) — 1 using Eq. (7) and found 
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TABLE 1 
M O 1 3 5 10 oo 
© 1.16 1.20 1.25 1.39 1.50 Loe 


that for a high Mach’s number, the velocity profile is 
approximately linear. The senior author,’ using this 
linear velocity profile, the integral relation between the 
friction and the momentum, and Eq. (7) found that 


_ Frictional force per unit width of plate 
(poU?/2) X length of plate 


1 
eq— 2 
1+ 5 M 


The dimensionless quantity 6 shown in Table 1 is a 
function of Mach’s number only. 
However, if Eq. (7a) is used, then 


It is evident that this linear approximation is not 

satisfactory for small values of Mach’s number. For 

M = 0, the case is the same as the Blasius solution’ 
for incompressible fluids for which 0 is 1.328. 

To solve the problem more rigorously, one has to 

resort to Eqs. (1) and (2). By introducing the stream 

function ¥ which is defined by 


p oy p 
Po ov Po 


Cy 


(8) 


Ox 


the equation of continuity, Eq. (2), is satisfied auto- 
matically. Now, if in Eq. (1) y is introduced as the 
independent variable as was done by von Mises* in 
his simplification of the boundary layer equation for 
incompressible fluids, and all terms are reduced to non- 
dimensional form then 


ou* 
where 
u* = u/U 
n* = n/L 
= (¥/UL)V pUL/us 
p* = p/po 
= (9a) 


and L is a convenient length, say length of the plate. 
Eq. (9) can be further simplified by introducing a 
new independent variable ¢ = y*/1/n*, 


then 


d * du* 
— 1 
a (« (10) 
This can be solved by the method of successive ap 
proximations. As p* and yu* are functions of tempera- 
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Fic. 2. Skin friction coefficients. (a) No heat trans- 
ferred to wall. (b) Wall temperature 1/, of free stream 
temperature. (c) von Karman’s first approximation. 


ture only as shown in Eqs. (6) and (7) or (7a) and the 
temperature is a function of u* then by starting with 
the known Blasius’ solution® the right-hand side of 
Eq. (10) can be expressed in terms of ¢. Therefore, 
one can write 


u*p*y* = 
Consequently, the solution of Eq. (10) is 


u* = cf Pat (11) 
o f 
where 
and C is determined by the boundary condition, 
1 
& lla 
C I 


A second approximation can be made based upon the 
value of u* obtained from Eq. (11). It has been found 
in the cases investigated that the third or fourth ap- 
proximation gives sufficient accuracy. 

Having computed the final ~*, the y corresponding 
to u* can be calculated from 


¥ VUpo/wox = dt/p*u* (12) 
Also the skin friction can be computed by the mo- 
mentum law, 
F 
* 2 (1 — u*)dg (13) 
9 VR 


Cy 


The velocity profile, the temperature distribution, 
and the frictional drag coefficient are calculated for 
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Fic. 3. 


Veiocity and temperature distribution when no 
heat is transferred to wall. 


different values of the Mach’s number of the free stream, 
for the case [(x — 1)/2]M? = (T7,,/To) — 1 using the 
approximate viscosity relation of Eq. (7a). The results 
are shown in Figs. 2 and 3. The velocity profiles for 
high speeds are very nearly linear, but it can be seen 
that the wall temperature for greater Mach’s numbers 
is very high. If the free stream temperature is 40°F., 
then the wall temperature will be 1600°F., 3620°F., 
6540°F., and 10,170°F. for Mach’s numbers of 4, 6, 8, 
and 10, respectively. Therefore, there is no doubt that 
the law of viscosity as expressed by Eq. (7a) will not hold. 
Also at such high temperatures, the heat transfer due 
to radiation cannot be neglected. Therefore, the re- 
sults for extreme Mach’s numbers are qualitative only. 

The change in the constant C;/R is appreciable, 
but not great. It decreases from 1.328 for M = 0 
to 0.975 for MZ = 10, or about 30 percent. However, 
for 0 <M <3 the change of the constant is very small. 

Fig. 2 also shows that Eq. (8a) which was obtained 
by using the linear approximation is fairly accurate for 
very high Mach’s numbers. 

As examples, consider first a projectile and second, a 
wingless sounding rocket. Taking the diameter of the 
projectile to be 6 in., the length 24 in., the velocity 
1500 ft./sec. and the altitude 32,800 ft. (10 km.), then 
the Reynolds Number based on the total length is 
7.86 X 10° and the speed ratio is 1.52. From Fig. 2 
the skin friction coefficient is 


Cy = (1.286 X 107*)/V7.86 = 0.000459. 


Changing the skin friction coefficient (based on the 
skin area) to the drag coefficient (based on the maxi- 
mum cross-section), one obtains 


Coy = 0.0055. 
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Fic. 4. Velocity and temperature distribution when the 
wall temperature is '/; of the free stream temperature. 


The drag coefficient due to wave formation takén from 
Kent’s experiments" is 


Cp,, = 0.190. 


Therefore, the ratio of skin friction to wave resistance 
is 0.0055/0.190 = 0.029. 

However, the ratio is greatly changed in the case of 
the rocket. Taking the diameter of the rocket to be 9 
in., the length 8 ft., and the altitude of flight 50 km.* 
(164,000 ft.), the velocity 3400 ft./sec., then the Rey- 
nolds Number based on a density ratio at that altitude 
of 0.00067 and temperature 25°C. (deduced from data 
on meteors) is 6.14 X 10°, and the speed ratio is 3.00. 
From Fig. 2, the skin friction coefficient is 


Cy = (1.213 X 10-*)/V 11.4 = 0.00360. 
Then 
Co, = 0.123. 


* The hydrodynamic equation holds so long as the mean free 
path of the molecules is small in comparison with the thickness 
of the boundary layer. For this case the thickness of the 
boundary layer is zero at the nose, however, at a distance '/; of 
the length of the rocket it already amounts to 3.2 cm., while the 
calculated mean free path of the air molecules at the altitude 
considered is about 1.1 K 107~*cm. Hence it appears that even 
for this case the theory can be safely applied. This conclusion 
is substantiated by the experimental results of H. Ebert in 
“Darstellung der Strémungsvorginge von Gasen bei neidrigen 
Drucken mittels Reynoldsscher Zahlen,” Zeitschrift fiir Physik, Bd. 
85, S. 561-564, 1933. 


The drag coefficient due to wave formation from Kent's 
experiments’ is 


Cp,, = 0.100. 


Therefore, the ratio of skin friction and wave resistance 
is now 0.123/0.100 = 1.23. If the boundary layer js 
partly turbulent, the ratio will be even greater. This 
shows clearly the importance of skin friction in the 
case of a slender body moving with high speed in ex- 
tremely rarefied air. It also disproves the belief that 
wave resistance would always be the predominating 
part in the total drag of a body moving with a velocity 
higher than that of sound. The reason underlying this 
fact can be easily understood when one recalls that the 
wave resistance of a body is approximately directly 
proportional to the velocity, while the skin friction is 
proportional to the velocity raised to a power between 
1.5 and 2. Therefore, the ratio of skin friction to wave 
resistance increases with the speed. With very high 
velocities and high kinematic viscosity, the wave re- 
sistance may even be a negligible portion of the total 
drag of the body. 


II 


In order to point out the thermodynamic aspect of the 
problem two cases will be considered: the flow of a hot 
fluid along a surface which is kept at a constant tempera- 
ture inferior to that of the fluid, and the case of a hot wall 
cooled by a fluid of lower temperature. The problems 
treated in this part have been discussed before in two 
very interesting papers by L. Crocco.**® He especially 
gives an elegant treatment of the cooling problem in 
the case of very high velocities (‘‘Hyperaviation’’). 
The authors feel that their treatment is somewhat more 
general and extended than Crocco’s previous analysis. 

An interesting general relation between the heat 
transferred through the wall and the frictional drag can 
be obtained using the assumption that Prandtl’s num- 
ber, 1.e., the ratio c,u/d, is equal to unity. The same 
assumption was used also in the previous calculations. 
It is remarkable that the relation holds also as well for 
laminar as for turbulent flow. The heat flow g per unit 
time and unit area of the wall surface is 


q = Av(OT/Oy) 
and the frictional drag 7 per unit area is 
T = wy, (Ou/Ov),, 


Using Eq. (4) the ratio g/7 can be calculated from thie 
relation 


where 7) is the absolute temperature, and LU the velocity 


of the fluid in the free stream, T,,. the absolute tempera- 


ture at the wall, ,. and yu, are the heat conduction and 
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viscosity coefficients of the fluid corresponding to the 


‘ wall temperature and .\ denotes Mach’s number. 


Substituting .1/ = 0 one obtains from Eq. (14) 


q To Te _ — Te (15) 
U pwU 
This is the relation known as Prandtl’s or G. I. Taylor's 
formula, first discovered by O. Reynolds. Hence Eq. 
(14) gives the correction of this result for compressi- 
bility effects. 

In the case 7) > 7),, 7.e., when the wall is colder than 
the free stream, the effect of compressibility is to in- 
crease the heat transferred through the wall. However, 
it would be erroneous to interpret this result as an 
“improvement” in cooling because at high speed the 
heat produced in the boundary layer is of the same order 
as the heat transferred through the wall. In order to 
determine the efficiency of the cooling a complete 
heat balance must be made. For this purpose Eq. (14) 
does not give sufficient information and the velocity 
and the temperature distribution in the boundary laver 
must be computed. Such calculations were carried out 
for the particular assumption = 75/4, 1.e., for the 
particular case in which the absolute temperature of 
the wall is kept constant at a value equal to one-fourth 
of the temperature of the hot fluid. With the same as- 
sumption for the variation of w as in Part I, the re- 
sults shown in Fig. 2 and Fig. 4 were obtained. ‘The 
variation of C;./R with / is similar to that obtained 
in the case without heat conduction through the wall. 
Also the highest temperature in the boundary layer is 
very high for extreme Mach numbers. However, the 
temperature maximum occurs some distance from the 
wall. 

The heat transferred from the boundary layer to the 
wall can be calculated as follows: 

The initial slope of the velocity profile is equal to 


ov w L V/n* bw 4 


By differentiation of Eq. (4) the relation between the 
velocity slope and the temperature gradient can be ob- 
tained. Using Eq. (7a) and substituting Eq. (16) into 
Eq. (5) then 

= K(To/R/4Ly/n*) (17) 
where 


K = (4°-76/2)(0.75 + [(« — 1)/2]M)V/RC, 


Therefore, the heat transferred to a strip of unit width 
of the wall of length L per unit time is equal to 


or approximately 
Q Ky. To/R (18) 


where K is given in Table 2. 
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Heat balance when the wall temperature is '/, of 
the free stream temperature. 


Fic. 5. 


TABLE 2 
VW K 
0 
1.93 
2 3.12 
10.53 
10 33.98 


The total heat balance is shown in Fig. 5. The ‘‘dis- 
sipation”’ curve represents in dimensionless form the 
heat produced by friction per unit time and unit width 
of the plate. The lower curve shows the increase (or 
decrease) of the heat content per unit time and unit 
width. The difference of the ordinates corresponds to 
the heat transferred through the wall. It is seen that 
cooling takes place for ./ < 2.6. Beyond this limit more 
heat is produced by friction than the amount which 
can be transferred to the wall and, as a matter of fact, 
the fluid is heated. 

In the case 7, > 7), 7.e., when the wall is hotter 
than the free stream, the ratio between the heat trans- 
fer and the drag decreases with increasing Mach’s 
number. This is shown in Fig. 6 where the ordinate 
represents the ratio between g/r with compressibility 
effect (according to Eq. (14)) to g/r without com- 
pressibility effect (according to Eq. (15)). The calcu- 
lation was carried out for a gas temperature of — 55°F. 
and a wall temperature of 1S0°F. and 300°F. It is 
seen that there is no cooling in the former case for M = 
1.69 and in the latter case for 1 = 2.08. However, 
the decrease of cooling efficiency is appreciable even 
at much lower speeds. This emphasizes the benefit 


| 
| 
= 
is 40 
a 
le | VA 
- 
ig 
y 
Is 
e 
h 
| 
| 
; 
) 
| 
Ou U VR Mo R C; 
: (16) 
| 
| 
/ 
) $ 


232 JOURNAL OF 
10 
| -55°F 
> 


Reduced Heat Trans 
° 
> 


/0 20 


02 


04 06 0. 10 
Mm 


(4 


Fic. 6. The effect of high speed on cooling efficiency. 

of the reduction of the speed of cooling air and the 
relatively poor efficiency of cooling surfaces exposed 
directly to a high-speed airstream. The curves in Fig. 
6 being derived from Eq. (14) apply to laminar as well 
as to turbulent motion. 
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Book Review 


Metal Airplane Structures, by Fravius E. Loupy; The Nor- 
man W. Henley Publishing Company, New York, 1938; 455 
pages, ill., $5.00. 

Metal Aircraft Construction, by M. LANGLEY; Pitman Pub- 
lishing Corporation, New York and London, 1937; 364 pages, 


ill., $4.50. 


It is exceedingly interesting to set these books side by side. 
While one book is written in England and the other in the United 
States, both of them have the similar objective of presenting a 
general survey of metal airplane construction. Both books are 
excellent, comprehensive, and excellently illustrated. The 
main difference between them is that Mr. Langley has restricted 
himself mainly to illustration and discussion of various types 
of construction and shop processes, whereas Mr. Loudy presents 
less descriptive material, but discusses basic elements of struc- 
ture without, however, going into detail stress analysis. Mr. 
Langley’s book is to all intents and purposes an encyclopedia of 
metal construction with illustrations and descriptions of almost 
every well known metal machine built either in Europe or in the 
United States. 

Mr. Loudy’s book, first of all, includes a study of structural 
elements, such as tie rods, struts, channels, columns, composite 
sections, extruded shapes, etc. For all such sections data are 
given in regard to area, weight, etc. For struts various charts 
are given for various strength columns. Similar charts are given 
for characteristics of aluminum channels. Thus this first section 
of Mr. Loudy’s book provides good reference to common struc- 
tural elements. The chapter on Welded Joints presents infor- 
mation which is generally well known, but nevertheless is 
interesting to have in one volume. In discussing Riveted 
Joints, the author gives formulae for efficiency of riveted joints, 
takes care of eccentric calculations in loads, discusses bearing 
values, and presents some very interesting illustrations of riveted 
spars and of riveted joints. Special rivets such as the Thompson 
tubular rivet and the Iseman rivet receive attention. In dealing 
with Stressed Skin Design, Mr. Loudy gives a variety of charts 
showing maximum loads for plates and in general discusses the 
basic principles of stressed skin design in useful fashion. Here 
again data are collected from a variety of sources and useful ref- 
erences are provided. The same remarks apply to the treatment 


given to corrugated aluminum sheets. Due attention is given to 
tension field beams and the pertinent formulae are nicely collected 
Here again the principles are set forth in clear 
Thereafter follows the review of the 
strength tests on thin wall duraluminum cylinders. From Chap- 
ter 7 and onwards, the material becomes descriptive. Various 
types of wings are briefly discussed, a classification of metal wing 
beams and spars is made, tests of torsional characteristics, etc., 
are presented. The chapter on Fuselage Design is also mainly 
descriptive with summaries of tests of fuselage deflection reported. 
A fairly complete bibliography is appended. 

Mr. Loudy’s book adds very little to the present knowledge of 
the subject nor is it a treatise on stress analysis. Nevertheless, 
it will be decidedly useful for reference, and to provide the de- 
signer or constructor or student with a general survey of metal 


for reference. 
and concise fashion. 


construction. 

On the other hand Mr. Langley’s book makes no attempt at 
structural analysis nor at the presentation of test data. The 
strength of Mr. Langley’s book lies in the large number of dia- 
grams and illustrations throughout the book. Also, construc- 
tion aspects are given due weight. For example in describing a 
Gloster rib, full details are given for its shop handling and almost 
everywhere in the book description of a design is followed by at 
least a hint of the shop methods employed. The book confirms 
the fact that British practice makes a far greater use of strip 
steel, but is in other respects closely analogous to American 
construction. 

Mr. Langley also includes a chapter on flying boat hulls and 
seaplane floats—wood and metal. A chapter on Subsidiary 
Structures deals with fins, retractable engine mounts, retractable 
chassis, strut ends, control columns, tank structures, etc. Work 
shop processes and detail design are always interesting. Prac- 
tice here is decidedly different from American practice, but may 
serve as an inspiration. Remarks on the spinning of fairing rings 
are instructive. Various processes of riveting are described. 
The appended tables and specifications relate of course to British 
practice. 

Both books evidently deserve a place in the library of the 
airplane designer or constructor. 

ALEXANDER KLEMIN 
New York University 
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Special Maps for Use with Loop Direction Finders 


JAMES M. COBURN, JR., New York University 


Presented at the Instruments Session, Sixth Annual Meeting, I. Ae. S. 
January 24, 1938 


HE radio range system with which the airways of 

the United States are equipped, while classed as an 
aid to navigation, in reality have not been so used to 
any great extent in the exact sense of the word “Navi- 
gation.”” Rather they have been used as an aid to con- 
tact flying in that the airplane proceeds from landmark 
to landmark just as it does in contact flying. In either 
case the exact position of the airplane and the ground 
speed is secured only by awaiting passage over a known 
landmark, the only difference being that in the case of 
“radio contact flying’’ the landmark is aural rather than 
visual. Position at other times, when following a 
defined airway, may be determined with satisfactory 
accuracy by the use of dead reckoning. 

The addition of rotatable direction finding loops to 
the equipment of the airplane, whether they be of the 
so-called radio compass type or of the older and more 
simple aural direction finding loop, offers possibilities 
of a technique more nearly approaching the true sense of 
the word “Navigation” in that they provide a means 
of determining position at any time, whether on or off 
an established airway, without waiting for the moment 
of passage over a known landmark. 

However, great care must be exercised in the intro- 
duction of new methods to the navigational technique 
of flight personnel. The complexities of operation of 
the modern airliner, coupled with the exacting nature 
of airways traffic control regulations, place a burden 
upon flight personnel which permit of little addition. 
Of a consequence, any change in routine must introduce 
a minimum of complications if it is to receive whole 
hearted acceptance. Unwelcome reception by flight 
personnel may easily cause resistance to the acquiring 
of the necessary familiarity with the equipment and 
technique, which may be acquired only in good weather, 
to the end that confidence in results secured when 
really needed might be lacking. 

Of primary importance, then, in the introduction of 
navigational technique, is the development of equip- 
ments and methods which will impose a minimum of 
additional complications, either mental or mechanical, 
upon flight personnel, all operations to be performed 
by existing personnel without leaving their seats or 
diverting their attention more than momentarily from 
their other duties. 

With these objectives in mind, before investigating 
the possibilities of equipment and methods which 


will meet them, it might be well to survey the accepted 
method of securing a radio position, or fix, using a stand- 
ard direction finder of either radio compass or aural type 
and the well known Department of Commerce Aero- 
nautical Sectional Charts. This is the method which 
has become standardized from many years use at sea 
for securing positions from both radio and visual bear- 
ings. 

The successive steps are as follows: 

(1) Tune the loop receiver to the frequency of the 
station upon which the bearing is desired and wait for 
positive identification. It may be necessary to rotate 
the loop to a position where maximum signal strength 
is secured. 

(2) Rotate the loop to secure a null indication in the 
ear phones, or zero pointer reading on the left-right in- 
dicator. 

(3) The angular indication of the pointer of the 
bearing indicator dial is read to secure the apparent 
relative bearing. The relative bearing is the angular 
difference between the fore and aft axis of the plane 
and a line joining the position of the plane and the 
radio station. 

(4) From the calibration curve of the loop installa- 
tion determine the corrected relative bearing corre- 
sponding to the apparent bearing. 

(5) Read the magnetic compass to determine the 
compass heading of the airplane or the compass bearing 
of the fore and aft axis of the plane. 

(6) Convert compass heading to magnetic head- 
ing by adding easterly deviation or subtracting westerly. 
Compass correction card must be consulted for these 
values. 

(7) Convert magnetic heading to true heading by 
adding easterly variation or subtracting westerly. To 
determine the direction and amount of the variation, 
consult the isogonic lines on the sectional chart at the 
estimated position of the plane. 

(8) The corrected relative bearing is added to the 
true heading of the airplane to secure the true bearing 
of the radio station from the plane. 

(9) Add or subtract 180° to or from the true bear- 
ing to secure the reciprocal true bearing or the true 
bearing of the plane from the radio station. 

(10) Establish true north through the radio sta- 
tion on a map including the radio station and the 
estimated position of the plane. 
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(11) With a protractor, the station in the center 
and aligned on the true north, lay off the reciprocal 
bearing. 

(12) With a straight edge, draw line of position on 
the reciprocal bearing through the station. 

(13 to 24) Repeat using a second station. The 
intersection of the two lines of position is the approxi- 
mate position of the plane. 

(25 to 36) Repeat using a third station if in any 
doubt about results obtained from two. 


DISADVANTAGES OF MARINE METHODS APPLIED TO 
AIRLINERS 


The difficulties of securing accurate results by follow- 
ing the foregoing procedure in flight may be easily 
imagined. Even among seafarers of many years ex- 
perience cases continually arise where deviation and 
variation are applied in the wrong way, by addition 
where there should have been subtraction and vice 
versa. With many steps to the routine, one or more of 
them may be forgotten easily, particularly in times 
of mental strain and it is in these times that accurate 
results are most needed. 

The disadvantage which precludes the possibility of 
the every day use of this method of position finding, 
however, comes to light when it is attempted to plot 
reciprocal bearings on the Aeronautical Sectional 
Charts within the confines of the cockpit. To do 
accurate plotting on any map requires that it be 
spread out upon a flat surface. In the marine service 
a chart table is provided for this purpose. Aero- 
nautical Sectional Charts are of such a size that to 
spread out one of them would require a chart table of a 
size that is impossible in the already overcrowded cock- 
pit. When it is considered that two or more of these- 
maps might have to be joined together to cover the area 
that normally could be expected to be required in 
practice, the impossibility of the situation becomes 
evident. 


METHODS OF SIMPLIFICATION 


Simplification which will eliminate one step may be 
accomplished by incorporating the calibration of the 
loop in the bearing indicator dial. This means that a 
special dial must be made with a scale distorted so that 
the reading secured will be the corrected relative bear- 
ing (see Fig. 1). Ifa fleet of airplanes is to be equipped 
with direction finding loops mounted in exactly the same 
place and with identical receivers, calibration may be 
made in flight with one plane to secure the data for the 
constructing of the special dials with which to equip 
the whole fleet. This provides a very worthwhile 
simplification. 

Additional simplification which makes this method of 
securing positions of practical value to the transport 
pilot rests entirely with the maps to be used and it is 
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Loop bearing indicator with distorted dial incor- 
In this case relative bearings, instead 
of reading zero to 360°, read zero to 180° right or left. 


Fic. 1. 
porating calibration. 


believed that this is answered by maps which were 
devised by Peter H. Redpath, navigational expert, 
for this purpose (see Fig. 2). 

A scale of 1: 2,000,000 permits folding or cutting of a 
map into sections, easily handled by the pilot without 
leaving his seat. Sufficient area would be covered to 
include on one section the radio range stations on the 
airway to be used as well as those within useable range 
on adjacent airways on both sides. The size of this 
folded or cut sheet allows plotting without a chart table 
or special plotting board, the back of a navigation kit 
or a clip board for holding log sheets being all that is 
required. 

The Lambert Conformal projection was chosen for 
the reason that a straight line on the map closely repre- 
sents a great circle and it is furthermore the projection 
which is used in the Aeronautical Charts. 

Cities and larger towns, with their names, rivers, 
lakes, and other visual landmarks were incorporated to 
eliminate the necessity for consulting other maps for 
this information. Elevations of the highest peaks and 
contours of elevations above normal cruising altitudes 
were included for warning purposes. 

The location of radio range stations only, and not 
commercial broadcast stations is given since radio 
range stations identify themselves every few seconds 
whereas it may be necessary to wait as much as fifteen 
minutes for positive identification of a broadcast sta- 
tion. Consequently, under present conditions of in- 
frequent identification, commercial broadcast stations 
are of little value. 

A feature which greatly simplified the mechanics 
of plotting, by eliminating the necessity of using a pro- 
tractor, was the printing of a compass rose around each 
radio range station. Since this feature also made un- 
necessary the establishing of the true North through 
the station, two steps of the previously accepted routine 
were at once eliminated. 
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Another step was eliminated by orienting these 
compass roses on the magnetic North rather than the 
true North. With the roses thus oriented, the con- 
version of magnetic to true bearings is unnecessary. 
It may be argued that the use of the true North as a 
reference is necessary as the variation is continuously 
changing and normally will not, consequently, have the 
same value at the radio station and at the point from 
which the bearing is taken. This is quite true were it 
necessary or desirable to work to extreme accuracy. 
When accuracy is desired the correction for the differ- 
ence in variation may be secured from the isogonic 
lines on the map, together with the correction for the 
error occasioned by the convergence of the meridians. 
However, with the abundance of radio range stations, 
the distance from any two workable stations that 
might be selected will probably be such that the cor- 
rection is negligible. Nevertheless, the application of 
these corrections will be discussed later. 

To eliminate another step, that of securing the re- 
ciprocal bearing, the magnetic compass rose, instead 
of having its zero degree position placed on the mag- 
netic North was made a reciprocal rose by placing the 
180° position to the magnetic North. Thus the time 
necessary to add or subtract 180° to or from the 
magnetic bearing is saved and the attendant possibility 
of error is removed. 

With the simplifications introduced by these changes 
the routine of securing a radio fix appears very differ- 
ent, the steps being as follows: 


(1) Tune the radio to the desired station. 

(2) Rotate loop to secure a null indication. 

(3) Read indication of the pointer on the bearing 
indicator. 

(4) Add the relative bearing to the compass heading. 

(5) Add or subtract the deviation. 

(6) With a straight edge, draw line of position 


through the magnetic bearing on the reciprocal rose 
around the selected station. 
(7) Repeat on second station. 


The intersection of two lines of position is, of course, 
the approximate position of the airplane. In all but a 
few cases, where extreme accuracy is desired, or the ap- 
proximate position is found to be several hundred miles 
distant from one or both of the stations, this approxi- 
mate position will be sufficiently accurate for all prac- 
tical purposes. However, from this approximate posi- 
tion it is very easy to work out a corrected position in 
cases where it is necessary or desired, as will be seen 
later. 

Maps constructed along these lines were produced 
under the direction of Peter H. Redpath by Transcon- 
tinental & Western Air in the form of strips with the 
main TWA routes in the center. The scale, | :2,000,- 
000, is such that radio range stations on airways both 
North and South are included on a strip fourteen inches 
wide. 

Instructions were issued to flight personnel outlining 
the exact procedure to be followed in obtaining a fix 
using these maps, as a supplement to a navigation 
course covering technical aspects of navigation by ra- 
dio bearings in general and the various uses to which 
directional radio apparatus aboard a plane might be ap- 
plied. Flight personnel were instructed to plot fixes at 
least once on every trip, checking the accuracy by 
visual landmarks where possible, to report upon the 
degree of accuracy and the time required to secure a fix 
by this method. 


REPORTED RESULTS 


The reports of results obtained were most satisfac- 
tory in nearly all cases. Positions were determined 
with an error of from one and one-half to five miles 
within one minute when using nearby stations. Ac- 
curacy within these limits left little to be desired and it 
was at first thought that it would be unnecessary to 
burden flight personnel with rules for applying correc- 
tions for the difference in magnetic variation and for the 
convergence of the meridians. Conditions arose, how- 
ever, which made it advisable to include rules for these 
corrections in the instructions to flight personnel. 

With certain conditions of radio reception and again 
in certain localities, bearings were obtained upon dis- 
tant stations with greater ease than upon stations 
which were comparatively nearby. In these cases er- 
rors were reported of as much as 25 miles. Checking 
the results showed that this error was caused almost 
entirely by failure to allow for the difference in magnetic 
variation between the position of the airplane and the 
station upon which bearings were taken and the plotting 
error caused by the convergence of the meridians on the 
Lambert Conformal projection. 

Simple rules were, therefore, worked out which could 
easily be memorized or carried in convenient memoran- 
dum form. These rules are: 
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(1) When the approximate position of the airplane 
is West of the station, add six tenths of one degree for 
each degree difference in longitude and also add one de- 
gree for each degree difference in variation. 

(2) When the approximate position of the airplane 
is East of the station, subtract six tenths of one degree 
for each degree difference in longitude and also sub- 
tract one degree for each degree difference in variation. 


If, then, an airplane desires a routine position check 
from distant stations the foregoing rules are applied to 
the dead reckoning position, but if radio reception has 
been lacking for some time it may be a saving in time to 
plot the approximate position, determine the amount of 
the correction to apply from the relation of this ap- 
proximate position to the position of the stations upon 
which bearings were taken and then replot the actual 
position. 

As an example of the use of these maps, assume that 
an airplane is flying on a magnetic heading of 78° and 
is not following the beam course of an airway. The 
position of the airplane is not known and it is desired to 
secure a radio fix. 

Kansas City is tuned in at good volume so is selected 
as one of the triangulation points since the volume level 
indicates that it is not far distant. The loop is rotated 
for a null indication and the relative bearing is found to 
be 94°. 

Kirksville, Missouri, is next tuned in and the relative 
bearing is found to be 22°. 

Reducing these relative bearings to magnetic gives 
78° plus 94° or 172° on Kansas City and 78° plus 22° 
or 100° on Kirksville. 

Plotting these results, by use of the magnetic roses on 
the map, locates the airplane six miles southeast of 
Bedford, Iowa. (See Fig. 2.) 

After it was definitely determined that a satisfactory 
degree of accuracy could be obtained, recommendations 
were made to the Bureau of Air Commerce that maps 
of the entire country, constructed in this manner, be 
prepared so as to be available for use by pilots of planes 
equipped with directional radio apparatus. This pro- 
gram is now under way and Aeronautical Chart for 
Radio Direction Finding No. 24 D.F. is now available. 


ADDITIONAL REQUIREMENTS CONTRIBUTING TO 
ACCURATE RESULTS 


Up to this point, equipment and methods only have 
been studied and no attention has been paid to the hu- 
man element. Before short cuts may be safely em- 
ployed a thorough understanding of the underlying 
principles is necessary. Education, therefore, is es- 
sential in the theory of navigation by cross bearings, 
in the manipulation of the direction finding equipment 
itself, the taking of bearings, the solution of orientation 
problems and the many other uses to which directional 
radio apparatus may be put. 


The mere equipment of an airplane with a directional 
radio device and the proper maps in itself is no guaran- 
tee of additional safety. Rather a false sense of security 
may be engendered unless steps are taken to assure that 
flight personnel acquire this familiarity with the equip- 
ment and methods. 

As an instance, if the radio compass or aural direction 
finding loop is used as a simple homing device for fol- 
lowing an airway in the manner which is generally ad- 
vocated by the manufacturers of certain devices of this 
nature, there is an even chance that the plane will drift 
to the wrong side of the airway with attendant possi- 
bility of collision with traffic from the opposite direction. 
Manufacturer’s pamphlets, demonstrating in a dia- 
grammatic manner the curved flight path resulting from 
homing in a cross wind, always show the drift as being 
to the right of the center line of the airway. The fact 
that traffic in the reverse direction would be drifting to 
the left of the center line of the airway and in the way of 
traffic from the opposite direction, is strangely unno- 
ticed. 

This course of education should be standardized for 
all users of this equipment and they should be care- 
fully examined to determine their ability to navigate 
with this equipment alone, under all conditions that 
may arise, lest they become a distinct menace to them- 
selves and other traffic using the airways. 


ACCURACY REQUIREMENTS FOR AUXILIARY 
EQUIPMENT 


Just as important is the assurance that reliance may be 
placed on auxiliary equipment upon which the accuracy 
of radio bearings depends. Since it is primarily necessary 
to know the angular relationship of the longitudinal 
axis of the plane to the magnetic North, which can be 
secured from the magnetic compass only, and the rela- 
tive bearing of the radio station to the longitudinal 
axis, it is a requisite that the compensation of the mag- 
netic compass and the calibration of the loop installation 
be performed by thoroughly competent men who will 
certify them to be correct. 

With equipment and training such as has been out- 
lined, it has been found that position finding is not an 
operation steeped in mystery, to be attempted only in 
cases of dire necessity, but rather one that is welcomed 
by flight personnel as a navigational method to be regu- 
larly employed in daily operations and one which 
rapidly becomes indispensable. 

Obviously, if the addition of these special direction- 
finding maps coupled with simple training and educa- 
tion in the use of the combined equipment is all that is 
necessary to secure such results, these maps should be 
included, as fast as issued, in the equipment of all pilots 
using directional radio equipment and steps should be 
taken to assure a standard degree of proficiency in di- 
rectional radio navigation. 
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On the Stability Problem in Hydrodynamics, II 


C. L. PEKERIS, Massachusetts Institute of Technology 


(Received February 23, 1938) 


ABSTRACT 


T IS shown that the parabolic flow of a viscous 

fluid between two fixed walls is stable for small 
disturbances which have symmetrical streamlines. 

(1) The striking phenomenon of the breakdown of 
the laminar flow of a viscous fluid has so far received a 
theoretical explanation only in the case of fluid which 
is confined between rotating cylinders.! The in- 
stability which is found there is also known to exist 
for an inviscid fluid and the effect of viscosity appears 
to be to defer the breakdown until some critical value 
of Reynolds Number. The parabolic laminar flow be- 
tween fixed walls, which is the subject of the present 
discussion, was shown by Rayleigh to be stable in the 
limit of vanishing viscosity and small disturbances. 
It would thus appear that the instability which is ob- 
served is due to viscosity, a result which does not seem 
plausible. In fact, it is shown in this paper that 
in the particular case of symmetrical disturbances, 
stability persists for all Reynolds Numbers. The un- 
symmetrical disturbances have been discussed in a 
previous paper? for a limited range of small Reynolds 
Numbers, where it was shown that their damping is 
everywhere greater than for the symmetrical dis- 
turbances. It is thus not likely that the unsymmetrical 
disturbances will be unstable where the symmetrical 
disturbances are stable, although an explicit proof to 
that effect does not exist at present. If the above 
presumption is true, it must be concluded that the 
method of small oscillations is incapable of explaining 
the breakdown of the parabolic flow, the explanation 
for which must be looked for in the deviations of the 
actual conditions from those assumed in the theory. 
Possibly the breakdown can be brought about only by 
disturbances of finite instead of infinitesimally small 
amplitudes, as would seem to be indicated by the 
fact that in the ensuing turbulent flow the quadratic 
terms in the equations of motion play a significant 
role. Another possibility, which has been considered, 
is the inexistence of the parabolic flow at the intake 
and the possible instability of the boundary layer 
which is formed there. These factors involve even 
greater mathematical difficulties than the method of 
small oscillations and it would seem that their analysis 
will not be forthcoming until there is a satisfactory 
theory of turbulence. 

(2) Let the motion take place in the x-direction 
and be bounded by the planes y = +h. If the half- 


width of the channel / is taken as the unit of length, 
the undisturbed laminar velocity Up, which is assumed 


to be uniform along the x-direction is given by 


Uo 


(1) 


The unit of velocity is equal to the maximum velocity 
at the center. By the method of small oscillations it 
is assumed that a perturbing velocity-field is super- 
imposed on the above, which is of sufficiently small 
magnitude so that the squares of the velocities are 
negligibly small in comparison with their first powers. 
The above limitation is made not for physical reasons 
but rather for mathematical convenience, for other- 
wise the equations become non-linear and therefore 
untractable. Since the fluid is incompressible the 
velocities can be derived from a stream function 


u = O/dy,v = — Oy/Ox (2) 

where y is assumed to have the form 
= f(y) exp. ia(Ct — x) (3) 
y is thus periodic in x of wavelength \ = 27/a. C is 


in general complex = a + 7b, where a is the phase 
velocity and the motion is damped if 6 is positive, or 
unstable if } is negative. 

It is found? that f(y) satisfies the equation 


+ +i aR[2f + =0, (4) 


where the dots stand for differentiation with respect 
to y and R = Uh/», U being the absolute magnitude 
of the laminar velocity at the center and v the kine- 
matic viscosity. The condition of the vanishing of 
u and v at the walls implies that 


fQ) =f(-1) = f(-1) = 0 (5) 


Eq. (5) is satisfied only for particular values of C and 
the solutions of the system of Eqs. (4) and (5) for these 
values of C represent the possible oscillations of the 
parabolic flow which are periodic in x. From Eggs. 
(4) and (5) it is clear that f(y) is either even or odd, so 
that the perturbations have either symmetrical or anti- 
symmetrical streamlines. The stability problem re- 
duces then, in the above, to the finding of the proper- 
values C and especially to the determination of the 
signs of their imaginary parts b. If the C’s are all 
found to be positive the motion is proved to be stable 
for small disturbances, and unstable otherwise. 

It was shown in I that the b’s are all positive and the 
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TABLE 1 


Values of aRC = aRa + iaRb from the Third Approximation 


aR 10 100 1000 10000 100000 

a aRa aRb aRa aRb aRa aRb aRa aRbh aRa aRh 
0 3.52 9.48 | 2.16 21.65 0.21 21.60 0.021 21.60 0.0021 21.6 
0 4.74 27.74 | 54.63 19.98 | 555.48 20.40 | 5555.55 20.40 | 55555.6 20.40 
0 6.18 34.78 | 87.65 30.36 | 888.75 30.00 | 8888.87 30.00 | 88888.8 30.00 
0.5 3.86 9.16 | 9.61 21.48 71.90 21.40 716.69 21.40 7167 21.40 
0.5 | 4.94 27.41 | 55.11 19.56 | 562.63 20.03 | 5627.2 20.02 | 56273 20.02 
0.5 6.44 35.04 | 87.71 30.57 | 889.77 30.18 | 8899.2 30.19 | 88990.7 30.19 
0 4.71 8.68 | 24.84 22.17 | 216.35 21.89 | 2160.4 21.89 | 21589.6 21.89 
1.0 5.37 26.97 | 56.80 18.01 ] 585.99 18.80 | 5861.3 18.89 | 58613.1 18.81 
0 6.88 35.73 | 87.87 31.20 | 892.7 30.69 | 8929.41 30.60 | 89308.6 30.68 
0.4r | 5.20 8.60 | 31.92 23.44 | 283.90 22.84 | 2835.5 22.84 | 28346.1 22.78 
0.4r | 5.57 26.96 | 58.51 16.67 | 605.31 17.88 | 6054.6 17.90 | 60522 17.93 
0.4" | 7.07 36.95 | 87.99 31.68 | 895.05 31.06 | 8952.4 31.05 | 89557.7 31.08 
10 8.61 103 56.14 119.40 | 538.79 119.3 5385. 12 119.30 | 53850.2 119.31 
10 6.83 119 91.36 106.77 | 872.85 117.7 8715.83 118.08 | 87159.3 118.07 
10 8.49 134 91.79 129.20 | 981.36 118.4 9829 118.0 | 98290 118.0 

100 8.15 | 10018 55.5 10018 555 10018 5553.6 10018 55536 10018 

100 8.15 | 10018 91.4 10018 889.8 10018 8886.9 10018 88869 6 10018 

100 8.15 | 10018 97.5 10018 999 10018 9998 10018 99979 '10018 

motion therefore stable, for both the symmetrical and — 

Sife—fife = 0 (9) 


antisymmetrical types of disturbances for values of 
aR less than about 100. Similar results were ob- 
tained by a different method for the symmetrical 
disturbances by S. Goldstein.* It was further proved 
that stability exists for values of aR which are less than 
the smaller of the quantities 2a? and a‘. An inspection 
of the above theorem will show that it can be amplified 
to prove stability for values of aR which are equal to 
the smaller of the quantities 2a? and a‘. Thus, by 
an application of the above theorem stability can be 
ascertained for the range of values of a and aR which 
are to the left of the thick lines in Tables 1 and 2. 

It shall now be shown that for very long wavelengths 
such that a —> 0, but aR remains finite, the para- 
bolic flow is stable for symmetrical disturbances. 
Eq. (4) reduces then to 


+ iaR2f + (1-y?—O)f] 
which integrates to 
f + iaR[2yf + (1 — y? — Of] = (7) 


If f is even, the terms on the left of Eq. (7) are odd 
and therefore A = 0. One solution of Eq. (7) is 


fiw (8) 


0 (6) 


If fe is the second even solution of Eq. (7) one shall 
have a general even solution of the form Af; + Bfo. 
The boundary conditions, Eq. (5), imply that at y = 1 


But the left-hand side of Eq. (9) equals the value of the 
quantity in brackets in Eq. (7) at y = 1, if in the 
latter fo is written for f. Hence it follows that 


fll) =0 (10) 
Now multiply Eq. (6) by f, where f is the complex 
conjugate of f, and integrate from y = —l toy = +1. 
From the identities 


the boundary conditions and Eq. (10) one obtains 


ffdy]=0, (11) 


which shows that the imaginary value of C is positive 
for all values of aR. The range of stability has thus 
been extended to the values at the tops of Tables | 
and 2. 

(3) Methods of extending the stability investiga- 
tions to larger Reynolds Numbers were suggested in 
I, but these proved to be laborious and the more so 
the larger the values of R. It appeared, therefore, 
that a numerical method whose applicability is not 
limited to any definite range of values of R would be 
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TABLE 2 
Values of aRC = aRa + iaRb from the Fourth Approximation 
aR 10 100 1000 10000 100000 

a aRa aRb aRa aRb | a@Ra aRb aRa aRb aRa aR) 
0 3.80 9.67 11.63 37.92 0.956 36.58 0.0955 36.57 0.00956 36.57 
0 5.42 32.02 38.58 32.13 437 .23 33.82 4374.97 33.83 43749 .99 33.83 
0 5.63 55.40 72.18 32.84 750.07 36.20 7500.01 36. 27 75000 36.27 
0 6.41 62.91 90.1 57.1 936.75 53.40 9374.92 53.33 93750 53.33 
0.5 4.15 9.38 20.37 41.18 68.96 37.37 678.93 37.36 6788 . 24 37.37 
0.5 5.60 31.87 37.59 28.15 444.77 32.27 4451.30 32.29 44513 .30 32.27 
0.5 5.62 55.04 72.25 32.96 752.29 36.18 7522.36 36.24 | 75223.56 36.24 
0.5 6.66 63.10 90.1 57.1 937 .05 53.58 9378.17 53.51 93782.48 53.51 
1.0 4.74 8.99 29.26 48.45 200 . 82 41.52 1995.67 41.46 19955.48 41.46 
1.0 4.57 31.61 45.33 18.91 472.39 27 .26 4728.18 27 .33 47282 .25 27.34 
1.0 7.32 54.57 72.38 33.31 759.20 35.98 7591.99 36.04 75919 .85 36.04 
1.0 7.07 63.71 90.06 58.20 937 .96 54.11 9387 . 83 54.03 93879. 10 54.02 
0.4m | 5.40 8.89 | 32.97 50.83 | 256.90 44.92 | 2557.60 44.83 | 25574.69 39.71 
0.44 6.28 31.55 50.37 16.92 497 .49 24.14 4977 .06 24.25 | 49770.99 29.38 
0.44 5.65 54.62 72.38 33.31 764.76 35.69 7648.13 35.76 76481 .37 35.76 
7.25 64.20 90.06 58.2 938 .65 54.51 9395.20 54.41 93952 .78 54.41 
10 8.64 102.86 | 51.05 137.3 422.90 135.77 | 4221.16 135.75 | 42210.76 135.74 
10 6.90 121.44 | 71.76 107.1 734.98 131.80 | 7343.79 131.86 73437 .30 131.86 
10 6.63 145.89 | 91.35 130.6 975.32 134.65 | 9844.20 133.95 | 98449.97 133.95 
10 8.45 161.75 92.09 156.9 929 33 129.72 9216.16 130.38 92155 .02 130.38 

| | 

100 7.81 10032 59.88 10025 426.95 10033 | 4272.85 10009 42728.5 10009 

100 7.81 10032 78.11 10032 781.05 10032 7810.5 | 10032 | 78105 10032 

100 7.81 10032 81.54 10053 987 . 20 10045 9868 .68 | 10069 98686 .8 10069 

100 7.81 | 10032 92.9 | 10018 929 10018 9290 10018 |: 92900 |10018 


preferable. Such is the method of central differences 
which was used successfully by Jeffreys* in a problem 
of thermal stability. It consists in subdividing the 
range of integration into a finite number of equal 
intervals and using the finite differences for the differ- 
entials in Eq. (6). Let the interval from —1 to +1 
be divided into 2” equal parts so that the length of 
one interval d equal 1/n. Denote the values of f at 
some given point by Zp, and at the ends of the neigh- 
boring intervals by the scheme 


One then has for the first, second, and fourth differ- 
entials of f the expressions 
Z= (Z; Z= a? (Z; — 2Z) + Z-1), 
] 
= a — + — 4Z-1 + Z-2) (12) 


The procedure will be illustrated for the case n = 2, 
d = 1/2. Denote the values of f at y = 0, and 1/2, 
respectively, by a and b. The value of f at y = 1 


must be zero by Eq. (5), and in order to make f(1) 
vanish, a value 0 is assigned to f(3/2) even though y = 
4/2 is outside the range of integration. 


Since an even 


solution is being considered, one has to assign for f 


at y = —1/2, — 1,—3/2, the values b, 6. Then at 
y =0 

f =a,f = 4(2b — 2a), f'Y = 16(— 8b + 6a) 
and at y = 1/2 


f =6,f = — 2b), = 16(8b — 4a). 
It follows now from Eq. (4) that 


16(— 8b + 6a) — Sa?(2b — 2a) + ata + 


taR[2a + (1 — C)(8b — Sa — a’a)] = 0, (13) 
16(Sb — 4a) — Sa?(a — 2b) + ath + 
ia R[2b + (*/; — C)(4a — 8b — a%)] = 0. (14) 


In order that Eqs. (13) and (14) be consistent, the 
determinant of the coefficients of a@ and ) must be 
made to vanish. 

Hence: 


{96 + 16a? + at + iaR[2 — (1 — C)(S + a®)]} X 
{128 + 16a? + at + iaR[2 — (°/, — C)(S + a)]} = 
{— 128 — l6a? + iaR XK S(1 — C)} X 


{— 64 — Sa? + iaR(3 — 4C)}. (15) 
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This is a second order equation in C which upon 
development and denoting aRC by x becomes 


x2(32 + 16a? + at) — x[aR(24 + 24a? + 1.75a*) + 
i(768 + 3520? + + 2a%)] — (4096 + 153602 + 
352a4 + + a’) + + 0.75a4) + 
iaR(256 + 2480? + + 1.7508) = 0. (16) 


The two roots of Eq. (16) which are complex and are 
functions of a and aR correspond, respectively, to the 
first mode of oscillation with no nodal lines and to the 
second mode of oscillation having one nodal line. 
These roots turn out to have positive imaginary parts 
for all values of the arguments a and aR. One can in 
fact dispense with the explicit solution for the roots and 
determine the sign of their imaginary parts directly 
from Eq. (16) by a well known algebraic operation. 
But the latter procedure, although offering a great 
time-saving in the subsequent approximations, does 
not reveal the degree of convergence of the successive 
approximations. 

Clearly if intervals of '/; and '/, are next used, 
one shall obtain instead of Eq. (16) equations of the 
third and fourth degree, respectively. The solution 
of these latter equations, which turned out to be ex- 
tremely laborious, is given in Tables 1 and 2. The 
degree of convergence of the solutions for the range of 
values of a and aR which are to the right of the thick 
lines are of primary interest. There, the first roots, 
for which the convergence should be best, differ in 
absolute magnitude by less than 10 percent. The 
imaginary values differ considerably, due to the fact 
that the real parts are so much larger and conse- 
quently a small error in the latter involves a larger 
error in the former. The first roots in the column for 
aR = 10 are in close agreement with the values ob- 
tained from the formulae in I. For a = 100 the 
imaginary parts are approximated to less than 1 per- 
cent already in first approximation, where d = 1. 
It did not seem worthwhile to proceed with a fifth 
approximation, especially since the evidence as to the 
sign of the imaginary parts is sufficiently convincing. 

(4) It is of interest to compare the oscillations 
of a viscous fluid with those of an inviscid fluid. This 
can be done in the limiting case of long wavelengths. 
The oscillations of an inviscid fluid are then governed 
by the equation 


F+(1—y — Cf =0 (17) 


subject to the condition that f(1) = 0. The solu- 
tion of Eq. (17) being (1 — y? — C), it follows that C = 


0. In the case of a viscous fluid it can be shown that, 
if a vanishes in such a way that aR —> 0 as R—,» 
©, the oscillations are damped as if the fluid were still, 
C then vanishes for all modes of oscillation like (a/R), 
thus approaching in the limit the value for an inviscid 
fluid. When a@ vanishes in such a way that aR —» 
© as R—+» o the values of the first roots at the tops 
of Tables 1 and 2 show a similar behavior. The other 
roots seem to remain finite, but the real parts are con- 
tinuously decreasing with the successive approxi- 
mations, so that it is not impossible, though unlikely, 
that they too approach zero. 

If absolute units of length and time are reverted to, 
y is of the form 


(CU 
f(y) exp. ta 


Now it seems from Tables 1 and 2 that with the ex- 
ception of the first root at a = 0, C approaches with 
increasing R a value (a + ib/R). 

Hence: 


aUt ibty «x 
h h 


(18) 

Now consider a disturbance of wavelength 27/a 
which is generated at the intake. According to Eq. 
(18) the time T it will take this wave to cross a given 
length L of the channel is of the order of L/aU. Dur- 
ing that time its amplitude will have decreased to 
exp. (—abvL/aU) = exp. (—abL/hR). Hence, as the 
Reynolds Number is increased the fraction by which the 
amplitude of a given small disturbance will be damped 
after crossing a given distance will decrease continuously. 
Presumably this tendency for increasing persistence 
with increasing R may turn to instability as the 
amplitude of the disturbance is increased. 

It would be of interest to apply the method of central 
differences to velocity profiles with inflection points, 
which from the work of Tollmien® are known to be 
unstable for an inviscid fluid. It should be expected 
that viscosity will affect stability for small Reynolds 
Numbers and that beyond a certain critical Reynolds 
Number the flow will become unstable. 
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On the Triple Origin of Air Forces* 
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IR forces have their origin partly in the shape of 
an obstacle obstructing the free motion of air, and 
partly in the condition at its surface; both causes gener- 
ally contributing together at the same time. A third 
and important item is the rate of change of these condi- 
tions giving rise to the creation of the air forces, so that 
the acceleration or change of flow pattern must be taken 
separately into account. 

That is well known, but not enough attention has 
sometimes been given to the fact that these airforce 
portions or components are of different dignity. The 
shape forces are of higher dignity and may be called 
noble, the surface forces are of less dignity and must be 
considered as base. This distinction is even more pro- 
nounced under accelerated conditions. Noble and base 
are to be taken in the sense used in old times by the al- 
chemists and still used by the modern chemists with 
respect to metals; gold being designated as a noble me- 
tal, and lead as a base metal. Modern chemistry also 
speaks of noble gases. Nobility and baseness are traits 
of human character reflected through all nature, which 
essential unity between humanity and the rest of nature 
is sometimes lost sight of in modern times. Formerly, 
people realized more clearly that inanimate nature de- 
picts and manifests the traits of human character. Be 
that as it may, there is no doubt but that the aerody- 
namic shape forces are much more positive, simple, pre- 
dictable and require less correction, in comparison with 
idealized theory, than surface forces. These can only 
be conditionally predicted, and require correction, and 
a less uniform and certain one, in applying theory. The 
acceleration shape forces in particular are almost cer- 
tain to make their appearance right away, nothing short 
of an actual compression of the air can retard them. 
The acceleration forces of surface effects, in contrast 
thereto, can hardly be predicted at all, and as far as 
they can be, are slow to develop and uncertain to realize 
full strength. 

The digest of the total air force, as measured in the 
laboratory or elsewhere by the designer, requires there- 
fore as an intermediate step its analytical separation 
into noble and base constituents, the same to be united 
again for the final engineering use of the information. 


*In view of the author’s preeminent place in the history of 
aeronautical science, his philosophical views, even though 
expressed in a somewhat unorthodox fashion, should be of general 
interest. The Editor. 


In that way, and in that way only, can one wind-tunnel 
reading be made to answer a thousand questions, and 
the wind-tunnel tests have to answer that many ques- 
tions lest wind tunnelling degenerate into fruitless 
weighing and air blowing. The wind-tunnel model has 
to be dimensioned previous to the test, not afterward. 
The question, incarnated in the model, and put to the 
tunnel, has itself to be sensible, if a useful and sensible 
answer is to be expected. It is physically impossible 
to test all combinations imaginable. Very often, the 
design has to be based on past tests and yet combines 
features never combined before in any test. Therefore, 
let wind-tunnel results be available in a form to be used 
broadly, quickly, easily, and still reliably. 

These general remarks may be illustrated by the dis- 
cussion of certain aspects of airfoil forces. Even in two- 
dimensional flow, the wing forces are decidedly of mixed 
origin. The profile or friction drag is beyond doubt a 
surface effect, and therefore of decidedly base nature. 
It can hardly be predicted by pure reasoning. It fails 
conspicuously to follow the legitimate square law of 
Newton; it rather flirts and plays with the less legiti- 
mate law of Reynolds, without even showing particular 
faithfulness to that. The profile drag permits itself to 
be lured off the straight path by every secondary and 
accidental influence. The physical conditions of the 
airfoil surface, the interference of adjacent objects, 
local disturbances in the air motion, and other influences 
sway this base and unpredictable item. 

With the pure lift of the two-dimensional flow, as 
contrasted to the lift of aerodynamic induction, things 
are already a little better. Kutta’s principle of a 
smooth trailing edge flow is indeed a surface effect, but 
only a small portion of the surface, the trailing end, is 
involved. The forward portion of the airfoil acts by its 
shape rather than by its surface to bring about the par- 
ticular condition at the trailing edge. The lift is there- 
fore less base than the profile drag, although it has not 
yet attained to nobility of the highest order. It is at 
least conditionally possible to precompute the lift by 
pure reasoning, and the correction needed to compute 
the actual lift from the theoretical one is not extreme. 
The zero lift direction of a wing section is even of higher 
nobility, and there is practical agreement between 
theory and practice in that respect. This is so because 
the absence of lift is predicated on the absence of sur- 
face effects calling for lift. As with humans, so also 
with surface effects: nonfeasance, doing nothing, sitting 
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still, can be more safely expected under conditions fa- 
vorable thereto, than positive action. 

Under simplest assumptions, the slope of the two- 
dimensional lift curve is predicted to be 27. During 
the war days, wind-tunnel tests prompted the author to 
substitute one tenth per degree for this theoretical value, 
a rule which has recently been confirmed and could not 
be improved by the most painstaking tests. One 
tenth corresponds to less than 27, so that there seems to 
be indicated a general lift deficiency of 10 to 11 per- 
cent. The theoretical lift fails to materialize at least 
by that much. Further progress of practical and useful 
aerodynamics seems to be tied up with looking a little 
closer into this question. Propeller blades seem to have 
a much steeper lift slope than ordinary airfoils. The 
actual deficiency of airfoils is often much larger than 
11 percent, because the theoretical slope of 27 holds 
only for infinitely thin sections. With thick sections, 
the above factor has to be multiplied by the sum of 
chord and thickness, in first approximation, rather than 
by the chord only. Thus, wing sections 10 percent 
thick would have a theoretical lift slope of 27 times 1.1. 
Its actual lift deficiency is therefore near to 20 percent 
rather than 10 percent. Suppose now that the lift de- 
ficiency is caused chiefly by the thickness of the tail 
portion, or behind parts, of the wing section, this being 
the place of base surface effects, rather than by the over- 
all thickness of the section. What would then be ex- 
pected from a thin section having a thick tail, so that 
the theoretical lift increase caused by overall thickness 
is absent, but the lift deficiency caused by tail thickness 
is present? A lift slope smaller than the almost gener- 
ally found one tenth per degree. This crocodile effect 
suggests itself easily to the mind weighing off baseness 
against nobility, in advance of any measuring. This 
very same thing has recently been verified experiment- 
ally The N.A.C.A. 0012-65 section is such a croco- 
dile section with a small body and a thick tail. Its 
lift slope is as much as 16 percent smaller than the al- 
most universally observed one tenth. Here then is a 
case where the secondary influences, which ordinarily 
do cancel, fail to cancel in consequence of unusual com- 
binations. This shows then that it pays to consider 
separately those items which ordinarily cancel each 
other, to separate what is separable, and to unite again 
at the proper time. This is a case where separation 
enables the advanced research worker to answer a ques- 
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tion about an unusual wing section from tests made 
with the usual wing sections. 

Applying the same logic to the question of airfoil 
equilibrium, let it be remembered that the moment 
about the point of application of the circulation lift is 
particularly noble, being in its entirety a shape effect. 
This central point of pure lift application is, with thin 
wing sections, at about 50 percent chord station. With 
thicker ones the central point tends to shift back, say 
by one fourth of the radius of curvature of the leading 
edge under ordinary conditions. The theoretical mo- 
ment about that central point is probably always 100 
percent realized. The lift is not that fully realized, as 
just discussed. Now, the distance of the important 
aerodynamic center from the central point is determined 
by the ratio of the moment slope to the lift slope. 
There results 25 percent in the simplest theoretical case, 
and hence the first choice of the aerodynamic center is 
50 percent minus 25 percent giving 25 percent. Ina 
representative case, the radius of curvature of the edge 
may be 4 percent, shifting the central point to 51 per- 
cent. The moment is not deficient, but smaller than 
the one corresponding to zero thickness by '/2 times 
0.1°, giving '/2 percent. Hence the aerodynamic cen- 
ter should be ahead of 51 percent by 25 percent X 
0.995/.9 or by 23.5 percent, giving about 23 percent. 
It actually is generally about at 24 percent. All these 
relations should ke carefully checked and the truth es- 
tablished. They are not final. 

Now suppose again that the airfoil has a crocodile 
section with a thick tail and slender body. The cen- 
tral point is then less displaced to the rear and closer to 
50 percent, but the distance of the aerodynamic cen- 
ter therefrom is increased, because the lift slope is de- 
creased. <A crocodile section should therefore have an 
unusually forward position of its aerodynamic center. 
The wind-tunnel test with the N.A.C.A. 0012-65 sec- 
tion again confirms this reasoning. Its aerodynamic 
center seems to be at 22 percent, closer to the leading 
edge than usual. 

The broad subject is far from being exhausted, but 
the space and time allotted to this paper is. These 
examples may be enough to illustrate the general prin- 
ciple. March divided and fight united is not only a 
good military rule of tactics, but it is also of highest value 
for obtaining fruits of scientific research which are 
really digestible by the engineer. 
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A Water Tank for Model Tests on the Motion of Airships in Gusts 
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SUMMARY 

A water tank and the auxiliary apparatus for determining the 
free motion of airships in gusts is described. The setup is in 
operation at the Guggenheim Airship Institute and a series of 
tests, financed by the Navy Bureau of Aeronautics, is being 
carried out. The method utilizes a water basin across a section 
of which a stream of water, simulating a gust, flows. An air- 
ship model immersed in the water and moved across the gust is 
free to turn and drift. The position and azimuth of the model 
are registered as functions of the time. The interpretation of 
the results in terms of forces aud moments acting on the ship 
are discussed. 


N AIRSHIP has a mean density equal to that of 
the air it displaces and so forms an integral part 
of the atmosphere in which it moves. Asa result it may 
be expected to carry out a considerable and complicated 
movement under the action of large scale gusts, 7.e., 
wind disturbances covering a volume comparable in size 
to that of the airship. This motion must be taken into 
account before an accurate determination of the forces 
acting on the ship can be obtained. 

This paper deals with a method for determining ex- 
perimentally the motion and, from this motion, the 
forces and pitching (yawing) moments acting while the 
airship flies through vertical (crosswind) atmospheric 
gusts. 

In order for the motion of a model through an arti- 
ficial gust to duplicate on model scale the motion of a 
full scale airship through an atmospheric gust, the fol- 
lowing conditions must be satisfied: (1) the model 
must have the same mean density as the fluid in which 
it moves, (2) the radius of gyration of the model about 
a transverse axis through its center of buoyancy must 
be / times that of the full scale airship, where / is the 
model scale, (3) the wave length or the width of the 
transition region from zero to full gust velocity of the 
artificial gust must be / times that of the atmospheric 
gust, and (4) the ratio of flight velocity to lateral veloc- 
ity in the gust must be the same for the model as for 
the full scale case. These conditions can be satisfied for 
a wide variety of gust structures by the method de- 
scribed below. 

As far as the author is aware this method represents 
the first experimental approach to the problem of de- 
termining gust forces on airships with a model which is 
free to turn and drift under the action of the gust. The 
investigation was suggested by the author and is being 
carried out at the Guggenheim Airship Institute with 
the financial support of the Navy Bureau of Aeronau- 


Fic. 1. Sketch of water tank. 


tics. This support and suggestions by members of the 
Bureau are gratefully acknowledged. 

The experiment is being carried out in water, since 
it would be extremely difficult to build a small scale 
model of an airship which would float in air and still be 
strong enough to withstand the gust forces. The method 
involves the use of a water basin (Fig. 1), across a 
section of which a lateral flow, simulating a gust, is 
maintained. The lateral flow is maintained in a closed 
channel intersecting the basin. The basin is 23 ft. 
long, eight ft. wide, and is filled with water to a depth of 
five ft. The channel in which the lateral flow is gen- 
erated is also eight ft. wide and five ft. deep. 

The airship model is suspended at its center of buoy- 
ancy from a carriage which is free to move laterally 
across the channel on tracks forming the cross members 
of a main carriage. The model is suspended at a depth 
of 2'/. ft. below the surface of the water. The main 
carriage and, with it, the model and cross carriage, are 
drawn along the basin in a direction perpendicular to the 
gust direction by means of weights acting over pulleys. 
The model is free to turn and drift under the action of 
the gust. The position of the center of buoyancy and 
the azimuth of the model are recorded by means of pens 
on a chart the movement of which is proportional to the 
movement of the model along the basin. Fig. 2 is a 
photograph of the setup. 

The model used is a 1/150 scale model of the AKRON- 
MACON, 62.8 in. in length, with a maximum diameter 
of 10.6in. It was cast in five parts of magnesium alloy 
with a wall thickness of about */s in. Weights were 
placed in the hollow space until the weight of the model 
plus cross carriage was equal to the weight of the water 
displaced by the model. The weights were distribu- 
ted inside the model so that the radius of gyration about 
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Cer. Fic. 2. Photograph of water tank showing carriage and 
airship model. 


a transverse axis through the center of buoyancy was 
1/150 that of the full scale airship. The moment of 
inertia of the AKRON-MACON was 73,000,000 kg. m. 
sec.” at sea level and 67,000,000 kg. m. sec.” at 3000 ft. 
Taking a mean value of 70,000,000 kg. m. sec.” the mo- 
ment of inertia of a 1/150 scale model which will float in 
water is 0.735 kg. m. sec.2» The moment of inertia of the 
model was determined by- comparing the period of oscil- 
lation of the model suspended by a wire at its center of 
buoyancy with that of a bar of known moment of in- 
ertia suspended by the same wire. 

Fig. 3 shows a gust velocity profile in a plane approxi- 
mately along the model path. The nose of the model 
strikes the gust along the line A, two feet from the wall 
of the basin. The model leaves the gust approximately 
along the line C, two feet from the far wall of the basin. 
In the intervening region the velocity is quite constant 
as is shown by observations along the lines B, B,, and 
B,. For this velocity profile the stern of the model is 
well into the gust before the nose starts to leave. The 
passage of the model into the gust represents in this case 
the passage of the airship through a wind shift such 
as is encountered particularly at cold fronts or during 
thunderstorms. 
¢ The main gust velocity for the profile shown is 2.5 in. 
; per sec. A model velocity of, say, 8.3 in. per sec. 
through the gust gives a ratio of gust velocity to model 
velocity of 0.3. The width of the transition region from 
zero to full gust velocity is about 7 in. The motion of 
the model under these conditions is, to scale, that of an 
airship flying through a windshift of 90°, the wind veloc- 
ity component normal to the initial flight path in- 
creasing from zero to 0.3 times the flight velocity over a 
distance of 150 X 7 in. or 87 ft. For a flight velocity 
of 100 ft. per sec. the corresponding gust velocity would 
be 30 ft. per sec. According to existing data '?* a gust 
of this magnitude with a transition width as low as 87 
ft. is near the upper limit of disturbances encountered 


in the free atmosphere. 
ee The Reynolds Numbers attained in these tests are 
necessarily low. Model velocities range from 6 to 18 
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Fic. 3. Velocity distribution across gust channel along 
approximate path of model. 


in. per sec., giving a range of Reynolds Numbers from 
60,000 to 160,000. 

The force in the gust direction and the moment of the 
forces about the center of buoyancy can be obtained 
from the cross and angular accelerations of the model, 
respectively. These accelerations can be found graphi- 
cally from the records of position and angle as a function 
of the time. 

One of the most important quantities derivable from 
these curves is the resultant quasi-steady angle of at- 
tack of the fins, particularly the maximum angle of at- 
tack attained for any particular value of the ratio of 
gust velocity to model velocity. This angle of attack 
with respect to the resultant flow follows directly from 
geometrical considerations taking into account the gust 
velocity and the motion of the model in determining the 
resultant flow. 

The severity of a gust of a given structure is deter- 
mined by the ratio of gust velocity to model velocity and 
by the width of the transition zone from zero to full gust 
velocity. In the set-up described both of these factors 
can be varied between wide limits. A variety of “gust 
structures’ can also be duplicated and by varying the 
width of the gust the effect of gust wave length on the 
forces experienced by the airship can be determined. 
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Determination of Profile Drag from Measurements in the Wake 
of a Body 


WILLIAM BOLLAY, Harvard University 


(Received March 2, 1938) 


INTRODUCTION 


ECENTLY the method of profile drag determina- 
tion by measurements in the wake of a body has 
received widespread interest by the aerodynamicists. 
It seems pertinent therefore at this time to give a simple 
deduction of the final working equations and to show 
clearly what the basic assumptions are upon which these 
equations are founded. 


FUNDAMENTAL PRINCIPLES 


When a real, viscous fluid flows past a body, the 
portion of it which comes in contact with the body is 
slowed down by friction. Thus there is a velocity 
deficiency in the wake of the body. Outside the wake 
the motion is not much affected. The problem is now 
to find a relation between this velocity deficiency and 
the drag which the body experiences. 

The total drag of a body is in general composed of 


two portions, the induced drag associated with a net . 


lift acting on the body and the profile drag due to the 
skin friction and separation losses. An attempt 
shall be made to determine expressions for these two 
drag components separately, making use of (1) the 
equation of continuity, (2) the momentum principle, 
and (3) Bernoulli’s equation along a streamline. 

Consider a control surface enclosing the body con- 
sisting of a cylinder with plane faces at the front and 
rear. Finally, let the radius of the cylinder tend to 
infinity. Then the continuity equation states that the 
net inflow into this cylinder is zero, 7.e. 


So Se 
where 


So plane face of cylindrical control surface ahead of 
body 

S; plane face of cylindrical control surface behind body 

S, cylindrical surface 

Uo, Uo, Wo are xX, y, 2 components of velocity at 
surface Sp 

U1, 01, W are x, y, components of velocity at surface S; 

u,, g, are the velocity components along and normal 
to cylindrical surface, respectively 

p is the air density 


Since the surfaces S; and S; are equal in area, the 
equation of continuity reduces to 


(1) 


Si 


The momentum principle states that the sum of 
the x-components of the forces acting on the fluid 
enclosed between the control surface and the body 
is equal to the outflow of momentum in the x-direction. 
This gives 


mi-as + f owa,-as + oui-as == 
So Sr 
So Sr Si 


where ~) and /; are the pressures acting at surfaces 
So and Sj, respectively, 7 is the shear force acting on 
the cylindrical surface and D is the total drag acting 


on the body. 
Since as r —> o, the velocity along the cylinder 


u,—> 


Sr S; 


In general, there is no shear force r acting on the 
cylindrical shell, and thus the equation reduces to 


D fv. + -dS - fv. + puj) — 
So Si 
—)-dS (2) 


or 
D = fv. — p)-dS + f — m)-dS (3) 
Si Si 


DEDUCTION OF EXPRESSION FOR INDUCED DRAG 


Eqs. (1) to (3) hold for an ideal fluid as well as 
for a real fluid. However, it is known that in an ideal 
fluid the total drag which can exist must be induced 
drag. Thus denoting with subscripts 7 the quantities 
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referring to the ideal fluid, the following relations are 
found. The equation of continuity becomes 


= f — u,)-dS = 0 (la) 
Sr Si 


This can be easily seen by considering the induced 
velocities due to a single horseshoe vortex and a doublet 
located within the control surface. Because of the 
symmetry there is no net influx or outflux over the 
infinite cylindrical boundary. The actual body can 
be made up by a superposition of doublets and horse- 
shoe vortex pairs, and thus Eq. (la) applies. 
The momentum principle gives 


D, = f + (2a) 
So Si 


or 
D, - fv. — py) dS fou —u,)'dS (3a) 
Sy Si 


The above integrals can be transformed by making 
use of the fact that Bernoulli's equation holds along a 
streamline even in vortex flow and thus the total 
head 


g = p + (p/2)(u? + + 


is the same along every streamline, namely, equal to 
that at Sp. Thus go = g,,. 
tions are rigorously correct. Now, however, the as- 
sumption is introduced that the plane Sp is far enough 
ahead of the body such that it experiences negligible 
induced velocities. Then hy = po + pu3/2 becomes 
equal to go, and one obtains 


D, = [go + (puj/2)]-dS — [es + (out,/2) — 
Se 
(0/2)(vi, + wi) ]-dS 
or 
D; = f + Wii) + f uss) 
(2b) 
or substituting in (3a), 
D; = f + dS — u,)?-dS 
Si S; 
(3b) 


In the limiting case where the surface S; is at an 
infinite distance behind the body which shall be de- 
noted by subscripts ©, “.; = mo (assuming straight 
trailing vortices) and 


D; = Po,) dS = [0/2003 + “dS 
Se So 
(4) 


Up to this point the deduc- , 


This equation shows that the induced drag is trans- 
mitted as a pressure difference between the front face 
and the rear face of the control surface. 


EXPRESSION FOR PROFILE DRAG 


The assumption made when considering the flow 
of a real fluid is that the induced drag is still given 
by the same formula as for an ideal fluid and that the 
remaining portion of the drag is profile drag. It is 
assumed that 7; = v,, and w, = w,; however, it re- 
mains to be checked by experiments whether p; = ,,. 

This leads to the following expressions for the profile 
drag, D, = D — D;; 


From Eqs. (2) and (2a): 


D, = — pi) -dS + u}) — 


Ss 
—u)-dS (5) 


D, = — p)-dS+ — uj) -dS — 
Si 


— m™)-dS (6) 


fw 21) -dS = f (p/2) (ui, + Vit wii) -dS— 
Si 


+ (0/2)(u2 + of + -dS= 
Si 


fi (p, — pi) -dS + (0/2)(ut, — -dS 
Si 


Everywhere outside the wake g, = g; and thus 


or using Eq. (la) 


Since go = g,; 


(go — gi) -dS 


fw — g)-dS = 
Wi 


where W, indicates the part of S, covered by the wake. 
Thus the above formulae transform into 


D, = J 2(g — g:) put — — 


Wy 
and 


D, = — g)-dS (ug — m)-dS — 
Wy Si 
(Pa — pi) dS (8) 


or 
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TABLE 1 
D, = (go — gi) + 

Wi x/c ACp, ACp, ACp, + ACp, 

(0/2) [(ui, — ui) — 2uo(u,, — m)]-dS (9) 

Ss, 1/4 — .00047 + .00013 — .00034 

In the case where the rear control surface is infinitely 1/2 — -00025 + .0001 — 00015 

1 — .00011 + .00007 — .00004 


far behind the body, Eqs. (3) and (4) give: 


D, = pu,,(uo — u,)-dS + J (boi — ba) (10) 
Se 


A formula analogous to Eq. (9) was first proposed 
by Betz! in 1925. For the second integral in Eq. 
(9) he carried out the integration only over W,, how- 
ever, and thus implicitly made the assumption that the 


integral 


T= (0/2) — ui) — 2uo(u,, — m)]-dS 
has a negligible value. Write 
uj, = + Au, + Aug 
Uy = uy + Au, + Aug + Au, 
where 
Au, = increase in u-velocity due to equivalent 
vortex 
Au, = increase in u-velocity due to equivalent 


doublet system representing the body. 

Au, = increase in u-velocity due to equivalent 
source which gives potential flow about a body 
having the form of the wake. 

Then 


I = | (p/2)- Au,:-(2Au, + Au, + 2Au,)-dS 


Si-Wy 
When the angle of inclination of the wing is small, the 


Au,: Au,-dS is approximately zero by 
S;-W, 
symmetry, since Au, is equal and opposite above and 
below the wake, while Au, has the same values. The 


integral fp - Au,: Au,-dS can be evaluated for any 

shape body by calculating for the corresponding doublet 
and source distributions. A rough estimate of the 
magnitude of this factor was made for a Rankine Oval. 
This gave the result that 


integral 


Jo: Au,: = ACp,: 


width of wake 
thickness of body 


to 
™ 


c = chord of body 
x = distance to section S, from trailing edge 


(p/2) - Au? - dS is calculated, 
S)-Wi 
substituting a single source at the center of the body, 
a second correction factor to the drag coefficient is 


found 


When similarly as above 


(29/c)? 
ACp2 
Sa (1/2) + (x/c) 
A numerical calculation for the case 2yo/c = 1/20, 


t/c = 1/8 gave the results shown in Table 1. 

These show that when measurements are made closer 
than about 1 chord-length behind the trailing edge 
of the body, the above correction factors should be 
applied.* The numerical values should probably be re- 
calculated where high accuracy is required for the 
particular airfoil shape used rather than using the 
above values derived on the basis of a Rankine Oval. 
Since an airfoil is not so full at the rear as the oval, 
somewhat smaller values would be expected. 

Jones? deduced a working formula’ for the profile 
drag starting his analysis with a modified form of Eq. 
(10), namely 


D, = fom. — 


i.e., he assumes p,; = p,. He transformed ihe 
above equation, making use of the fact that pu, -dS, 
= pu,-dS,; = dm is an element of mass which has 
suffered a loss in velocity (wu) — uw.) in unit time. 
Thus he writes 


D, =f —u,)-dS (11) 


* Since the submission of this paper a comment on the same 
subject by Muttray has appeared in the “Jahrbuch 1937 der 
Deutschen Luftfahrtforschung,’’ pp. 1-66. He states that the 
measured profile drags were found too high by about 12 percent 
at the trailing edge, by about 8 percent at a distance x/c=.10 
behind the trailing edge (compared to measurements a great 
distance behind the wing). In every case of a wing with finite 
thickness the Betz formula gave a variation of the profile drag 
with distance from the wing. Only in the case of a flat surface 
consisting of thin pieces of sandpaper glued back to back was 
the profile drag found to be independent of distance behind the 
wing. These results confirm the above conclusions very closely. 

t The same formula has also been deduced by Kramer and 
Doetsch, and had already previously been stated by Betz in 
1931. 
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where it is understood that (uw — u,) is associated 
with the particle which at infinity becomes pu,dS,. 
In order to determine u,, Jones makes the assump- 
tion that there is a hypothetical flow in which the total 
energy associated with a given particle remains con- 
stant in the wake in a similar fashion as it does for an 
ideal fluid with vortex flow. This is not entirely cor- 
rect since actually there is a small amount of dissipation 
of energy in the wake. The closer the plane of measure- 
ment S, is, however, to the trailing edge of the body, 
the less is likely to be the error. To get the final 


formula it is thus assumed that g, = g,, alsop, = 
Also for the 


Paw = Unt Vou and Ww, 
ideal fluid g..; = go and thus 


= + (0/2)[u, + + 
Bo = po + (pus/2) = poi + (0/2) + + 
Thus g = po + (pu’,/2) 
and V2/p V = 
u = V2/p — pi where hy = — p(v; + wi) 


Eq. (11) thus becomes 
D, = [2 V iy (Vg, Po V2, po) 
Si 


Since everywhere outside the wake g; = go this formula 


reduces to 
D, - V hy PLY V 2, Po) 
Wi (12) 


This formula thus has the great advantage over Betz’s 
that it should hold well when applied close to the trail- 
ing edge of the body without requiring any theoretical 
correction factors. It is to be pointed out that the 
surface W; should really be taken perpendicular to the 
free stream velocity and not perpendicular to the 
plane of the wing as Jones asserts. The two cases 
differ, however, only by a factor (1 — cos a)-D, 
where a is the angle of attack. 

The working formula corresponding to Eq. (7) or 
its equivalent forms Eqs. (8) or (9) becomes 


D, = — g)-dS— [2V g — — po — 
Vin — p)-dS+ 2D, (13) 


where AD, is a small correction factor which requires 
some theoretical calculations: 


AD, = (pis — Pi) + (p/2) - Au,- X 
WwW, 
(Au, + 2Au,)-dS 
Kramer and Doetsch® in their experiments at the 
D.V.L. simply assumed a formula of the type 


D, = (g — gi)-dS 
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They determined the mean value of the above integral 
by means of a whole array of total head tubes which 
were connected together so as to give the integrated 
value as a single reading. The constant K was deter- 
mined experimentally as a function of distance behind 
the wing. The above method requires still a consider- 
able amount of further investigation to prove that K is 
actually a universal function, independent of wing 
thickness and angle of attack. 

The main difficulties with the experimental deter- 
mination of the quantities required in working Eqs. 
(12) and (13) is first of all that the turbulence of the 
wake has a considerable effect especially upon the 
static pressure determination, as shown by the experi- 
ments of Schlichting. Secondly, in addition to the 
total head g, = pi + (p/2)(m? + v7 + wy”), the 
quantity h; = p, + (pu,?/2) is required as shown by 
Eqs. (12) and (13). The error caused by making the 
survey with a single total head tube to serve for both 
readings should be investigated, especially at the higher 
angles of attack when the quantity (o/2)(v? + w,?) 
assumes a larger value. Other sources of error such as 
the effective displacement of the center of a pitot tube 
in a pressure gradient have been investigated by Jones. 


SUMMARY AND RESULTS 


Formulae for the parasite drag from wake measure- 
ments corresponding to those by Betz and Jones have 
been deduced from basic principles. It is pointed out 
that Betz neglected quadratic terms of small velocities 
which, when measurements are carried out very close 
to the body may introduce appreciable errors. Jones's 
formula does not suffer from this defect; in fact it 
should hold best close to the body before the wake has 
suffered much energy dissipation. 

Both Jones and Betz neglected the quantity (o/2)- 
(v,? + w,”) in their formulae as small, which should be 
justified experimentally. Jones’s assumption of no 
energy loss in the wake appears reasonably good but 
requires also more thorough experimental verification. 
It should also be verified whether the pressure in the 
turbulent wake is the same as that calculated from a 
pure potential flow. It is pointed out that the plane 
of measurements in the wake should be normal to the 
free stream velocity, otherwise a small correction 
should be introduced. 
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Institute Notes 


NEw CORPORATE MEMBERS 


Five subsidiaries and the parent Bendix Aviation Corporation 
became Corporate Members of the Institute in March. The 
Lockheed Aircraft Corporation became a Corporate Member 


this month. 


TRIP OF PRESIDENT AND SECRETARY OF THE INSTITUTE 


The President and the Secretary of the Institute visited Chi- 
cago, St. Louis, Minneapolis, and Los Angeles during March and 
April. Major Gardner returned to New York through the cour- 
tesy of the Pan American Airways via Mexico City, Merida, 
and Havana, and flew from Miami to New York through the 
courtesy of Eastern Air Lines. A detailed account of their trip 
will appear in the next issue of the Journal. 


MEETINGS 


On March 8 about thirty members of the Institute gave a 
luncheon in honor of George E. Woods Humphery, Managing 
Director of the Imperial Airways, at the University Club, New 
York. 

On Wednesday, March 23, the Los Angeles Branch of the In- 
stitute, in cooperation with the American Society of Metals, 
presented a program of technical papers dealing with the use of 
metals in the aircraft industry. The first aeronautical session was 
held at the Biltmore Hotel in Los Angeles on Wednesday morning. 
The afternoon session was held in the Pan Pacific Auditorium. 
The program follows: 


Morning Session—9:30 a.m.—Biltmore Hotel 
Chairman 


C. L. Hrsert, Consolidated Aircraft Corporation 

J. R. Goipstern, Douglas Aircraft Company, ‘‘Metals in the 
Aircraft Industry.” 

ARTHUR W. WINSTON, Dow Chemical Company, ‘‘Recent De- 
velopments in the Application of Magnesium Alloys to Air- 
craft.” 

A. A. HANDLER AND L. W. Davis, Aluminum Company of 
America, ‘‘Design and Production of Aluminum Alloy Air- 
craft Forgings.” 

BRADLEY STOUGHTON, Lehigh University, ‘‘Engineering Alloy 
Steels (S. A. E. Specifications).” 


Afternoon Session—2:00 p.m.—Pan Pacific Auditorium 


Chairmen 


Howarp D. Houcuton, Douglas Aircraft Company 
E. E. SecuHiLER, California Institute of Technology 


WALDEMAR Navujoks, Steel Improvement & Forge Company, 
“Steel Forgings for Aircraft.” 

CaRL DE GANAHL, Fleetwings, Inc., ‘Stainless Steel for Air- 
craft.” 

Tom A. TrIpLett, Triplett and Barton Laboratories, ‘“X-Ray 
Inspection of Aeronautical Materials.” 


Plans are being completed by the Institute to hold national 
meetings in Los Angeles, Washington, and Chicago as well as in 
New York. Each meeting will have a special purpose and com- 
mand wide interest. Members and Branches should always 


remember in making plans that October 15 will always be cele- 
brated as the Anniversary of the founding of the Institute, and 
December 17 as the day of the Wright Brothers’ Lecture and Hon- 
ors Night. 


The Short Mayo composite aircraft. 


GIFTS TO THE INSTITUTE 


The donation of three book-cases and books by Major E. E. 
Aldrin is recorded on another page of the Journal with a list 
of the books he gave. Ernest W. Dichman also presented several 
hundred valuable pamphlets and reports. A. J. M. Hamon pre- 
sented F. T. Miller’s ‘‘The World in the Air’ and ‘‘Peru from the 
Air” by George L. Johnson. 

The Council acknowledges these donations with appreciation 
and hopes that other members will also wish to add to the rapidly 
growing library. 

Dr. Lyman J. Briggs, Director of the National Bureau of 
Standards, has presented a glass tube containing air from the 
stratosphere taken by Captain A. W. Stevens at an altitude of 
approximately 14 miles during the flight of Explorer II on Novem- 
ber 11, 1935. 

T. P. Wright, President of the Institute, presented a rib from 
the famous NC-4, the first flying boat to fly the Atlantic. 

G. E. Woods Humphery, Managing Director of the Imperial 
Airways, presented a large model of the composite aircraft 
designed by another Institute member, Major R. H. Mayo. 
Glenn L. Martin presented a large photographic mural of the 
China Clipper creating a whirlpool on one of her taxying tests. 

Charles Marcus gave a decorative stand holding a glass ball 
found by the Pan American staff at Midway. These balls are 
used by Japanese fisherman to float their nets. During storms 
the nets are blown to sea and sometimes these floats float across 
the Pacific and are found on the coast of California. 

The United Airlines presented a model of the DC-4. 

Another gift has been made to the Institute which will be an- 
nounced in the next issue of the Journal. 


AVIATION Morion PictuRE ANNUAL 


The Motion Picture Producers and Distributors of America, 
through the courtesy of Hon. Will H. Hays, are preparing for 
showing at Institute Meetings what will be known eventually 
as the Aviation Motion Picture Annual. All the aviation sub- 
jects that appear in the news films of six companies are sent each 
week to the Institute where they are joined and made into a 
continuous sound film. At the end of the year it will be edited 
and shown for the first time at the Seventh Annual Meeting of 
the Institute and later at Branch meetings. This courtesy on 
the part of the motion picture companies is in recognition of the 
great value aviation has been to the motion picture industry. 
This is but a preliminary announcement of the film of which 
members will hear more later in the year. 
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Fifth International Congress 
for Applied Mechanics, 1938 


The Fifth International Congress for Applied Mechanics will 
meet in Cambridge, Mass., September 12 to 16, 1938, at Harvard 
University and Massachusetts Institute of Technology. The 
Congress is a meeting of persons working in the field of Applied 
Mechanics who become members and may participate in its 
activities by payment of the registration fee of $5.00. A 
program of papers to be presented by members will be organized 
by selection from titles, with abstracts of 300 words, in the 
hands of the Secretariat by July 1, 1938. Such papers are lim- 
ited to 20 minutes for reading the paper and 10 minutes for dis- 
cussion. Papers presented at the Congress will be published 
in its Proceedings. 

General Lectures will be delivered by prominent scientists 
from England, France, Germany, Norway, and the United 
States on invitation of the Organizing Committee. 


The program will cover three main divisions of Applied 
Mechanics as follows: 

(1) Structures, elasticity, plasticity, fatigue, 
strength theory, crystal structure. 

(2) Hydro and aerodynamics, gasdynamics, hy- 
draulics, meteorology, water waves, heat transfer. 

(3) Dynamics of solids, vibration and sound, fric- 
tion and lubrication, wear and seizure. : 

Following the meeting at Cambridge, arrangements have 
been made to visit Washington (National Bureau of Standards), 
and Langley Field (National Advisory Committee for Aero- 
nautics). 

Dormitory facilities will be made available by Harvard 
University at a nominal charge. 


A bulletin giving further details will be mailed to those who 
address an inquiry to the 


Fifth International Congress for Applied Mechanics, 
Massachusetts Institute of Technology, 
Cambridge, Mass., U. S. A. 


Tu. von KARMAN 


ecretartes 
J. C. HunsaKer \s 
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MotTIOoN PIcTURE OF LANGLBY FIELD LABORATORIES 


The motion picture prepared for showing at meetings of the 
Institute by the National Advisory Committee for Aeronautics 
is now ready to be loaned to Branches. It is accompanied by 
comments prepared by the N.A.C.A. and takes about one hour 


to show. 


AERONAUTICAL REVIEWS 


Because of illness of the Editor of the Tecunicat Data D1- 
GEST, from which the Aeronautical Reviews published each month 
in the Journal are prepared, these reviews have not been avail- 
able for the April issue of the Journal. It is expected, however, 
that enlarged issues of the TecunicaL Data Dicest will be 
available for publication in the May issue of the Journal. 


Gift of Books to the Institute Library 


The following books have been presented to the Institute Library by Major E. E. Aldrin. 
The Council of the Institute thanks Major Aldrin for his generous gift of books that he has collected for many years. 


(continued from March 1938 issue) 


Precision of Measurements and Graphical Methods, by H. M. 
Goopwin; 1913. 

The Dynamics of Particles and of Rigid, Elastic and Fluid 
Bodies, 1912. Second Edition. 

Strength of Materials, by ARTHUR Mor.Ley; 1916. 
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History of Aeronautics, compiled by WiILL1AM B. GAMBLE; 
The New York Public Library, 1938; 325 pages, $2.00. 

This is the most useful bibliography of books, articles, and 
pamphlets dealing with the history of aeronautics that has been 
published. The author is Chief of the Science and Technology 
Division of the New York Public Library, and has spent four 
years compiling references to the aeronautical collection of over 
5000 volumes in the library. While the book gives references 
to outstanding events that occured after the Wright Brothers’ 
first flight in 1903, it is primarily concerned with every event be- 
fore that epoch making flight. 

There are many bibliographies on aeronautics but most of them 
presume that the enquirer knows the name of a book or its author. 
Mr. Gamble has arranged his book of 325 pages so that anyone 
wishing to secure information on any one of the many subjects 
that comprise aeronautical history may do so without much pre- 
vious knowledge of the subject. 

In addition to a careful division of subjects, short suggestive 
comment on almost every book or article gives the reader just 
the pertinent information he needs to select a reference. 

Even if the book had ended on page 254 it would have been of 
inestimable value, but the outstanding contribution of the author 
is the 71 pages of indexes. He gives complete references to all 
authors and subjects in the bibliography. But he goes much 
further. Under Bird Flight, for instance, a special bibliography 
is given and then the author lists all those experimenters who 
have written on this subject. Under Bombing Apparatus will be 
found a list of some of the earliest experiments in this field. 

It is extraordinary to find a librarian, who had no special in- 


terest in aviation, preparing such an exhaustive and helpful 
reference work as a part of his professional work. 

Every aeronautical writer, editor, and investigator will be 
greatly indebted to Mr. Gamble for making such a valuable 


contribution to aeronautical literature. 
LESTER D. GARDNER 


Jahrbuch der Deutschen Akademie der Luftfahrtforschung 
1937-1938; Deutsche Akademie der Luftfahrtforschung, Ber- 
lin, 1938; 321 pages. 

One of the most perfect examples of fine book making is this 
Yearbook of the Academy of Aeronautical Research containing 
biographies of all the distinguished members. The pictures of 
the members are reproduced in gravure and the biographies are an 
invaluable reference for securing information regarding many 
of the leading aeronautical scientists of the world. 

The Academy consists of only the most distinguished aero- 
nautical specialists in Germany and through careful selection 
includes as members a few leaders from many other countries. 
The care with which such selections have been made will be ap- 
parent from the list of American members which has been limited 
to Dr. Ames, Dr. Briggs, Dr. Durand, Dr. Hunsaker, Dr. Lewis, 
and Dr. Millikan. 

Another feature of the book is a list of all the papers that have 
been written by each member. The concluding chapters give 
the proceedings of the meeting of the Council held in April, 
1937. This book was presented to the Institute by Ministerial- 
rat A. Baeumker, who is Chancellor of the Academy, and it will 
make a valuable addition to the library of the Institute. 
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SUGGESTIONS FOR CONTRIBUTORS TO THE JOURNAL 


of the AERONAUTICAL SCIENCES 


The Institute of the Aeronautical Sciences invites 
both members and non-members from any country to 
submit papers for publication in the Journal of the 
Aeronautical Sciences. The Journal, following the 
practice of other scientific publications, does not pay for 
contributions. 


The following directions for the preparation of papers, 
if followed by authors, will save correspondence, avoid 
the return of papers for changes, minimize the work of 
preparation for the printer, and save the expense due 
to the charges made for “‘author’s corrections.” 


Manuscripts: Papers must be written in English, in original 
typewriting on one side only of white paper sheets, consecutively 
numbered, and be double or triple spaced with wide margins. 
Manuscripts should be prepared with great care so that they 
will be typographically accurate. Paragraphing should be 
given special attention. Papers should be written in the third 
person, reference to the writer being made as “the author.” 
Avoid the use of the words “I,” “‘we,” and “you.” Blueprint 
copies of papers are unacceptable as it is impossible to mark 
directions to the printer on them. Correcting, changing, or 
adding to papers after they are in type is costly. It is, there- 
fore, imperative that papers submitted be in final form. Typo- 
graphical errors may be corrected on proofs, but if authors wish 
to add material, they may do so at their own expense. In mail- 
ing, drawings may be rolled, but manuscripts should be sent 
flat. Send by first class mail (register if you wish for your 
own protection) to the Secretary, Institute of the Aeronautical 
Sciences, 5111 RCA Bldg., Rockefeller Center, New York City. 
All manuscripts will be examined by the Editorial Committee 
and by the Editor. Authors will be advised as promptly as 
possible (usually two to three weeks) whether the paper is ac- 
ceptable for publication. 


TrrLzs: The title of the paper should be brief. The name and 
initials of the author should be written as he prefers. The use 
of the full name of an author is advocated because of the fre- 
quent duplication of initials and surnames which sometimes 
makes it difficult to establish the identity of the author. This is 
particularly important for large annual indexing and abstract- 
ing services. All titles and degrees or honors are omitted. The 
name of the organization with which the author is associated 
should be placed after his name on the same line. The date on 
which the paper is received will be inserted by the Editor. 


SumMarRIES OR ABstTRACTsS: An abstract to be printed at the 
beginning should accompany each article. It should be ade- 
quate as an index and asa summary. It should contain a state- 
ment of major conclusions reached, since summaries in many 
cases constitute the only source of information used in compiling 
scientific reference indexes. Abstracts printed in other jour- 
nals, especially foreign, in most cases, consist of summaries 
from printed papers. The summary should explain as adequately 
as possible the major conclusions to a non-specialist in the 


subject. The summary should contain from 100 to 300 words, 


depending on the length of the paper. 


Sus-HeapINGs: Sub-headings should be inserted by the au- 
thor at frequent intervals. The work of editorial preparation 
will be simplified by the author providing many sub-headings. 
Owing to the breaking of columns and the insertion of illustra- 
tions, some of the sub-headings may have to be omitted. 


SHORTENING OF PaPERs: Some papers, at the end, fill in only a 
portion of a page. This leaves much wasted blank space as 
succeeding articles are started at the top of a page. Authors 
should indicate by notation on the left-hand side of the page 
what matter may be omitted when “run overs” occur. This 
request is important as the Journal cannot afford in the future, 
as it has in its earlier issues, to have blank half pages or more 
at the end of papers. 


MatrTer UsuaL_ty DELETED: Acknowledgments of assistance 
in preparation of paper, except by collaborators. Photographs 
or illustrations of little technical interest and not showing ad- 
vances in general practice. Too detailed tabular matter (gen- 
eral results of such tables may be included in the text). Lengthy 
descriptions of materials or processes or of preliminary experi- 
ments or theories which preceded final results; salient features 
only are of interest. 


REFERENCES AND FOOTNOTES: References should appear as 
footnotes only, numbered consecutively, grouped together at the 
end of the manuscript. The arrangement should be as follows; 
(for books) —! Durand, W. F., Aerodynamic Theory, Vol. 1, p. 23; 


Julius Springer, Berlin, 1934. (For magazines)— ‘Englund. 
C. R., Crawford, A. B., and Mumford, W. W., Some Results of a 
Study of Ultra-Shori-Wave Transmission Phenomenon, Proc. I. 
R. E., Vol. 20, No. 12, pp. 481 and 482, March, 1933. Please 
give Author, Title, Volume, page, publisher and date of publi- 
cation as indicated. Omission of one required fact causes much 
extra editorial work and possible inaccuracies. All references 
are grouped at the end of the article 


ILLUSTRATIONS: Illustrations should accompany manuscripts 
and each should always be referred to in the.text, preferably by 
number. Drawings or graphs should not be larger than 12 X 16 
inches, and must be made with jet black India Ink on white 
paper or tracing cloth, the latter being preferred. Do not use 
typewriter for lettering. The smallest lettering on 8 X 10 inch 
figures should be no less than '/, inch high. Cross-section 
paper (white with black lines) may be used, but should not have 
more than 4 lines per inch. If finer ruled paper is used, the 
major division lines should be drawn in with black ink, omit- 
ting the finer divisions. In the case of finely ruled paper only 
blue-lined paper can be accepted. Tracing paper and blue- 
prints are not acceptable. Lettering and all markings must 
be large enough to be readable after reduction. Mail rolled or 
flat, never fold. Drawings which cannot be reproduced (includ- 
ing pencil drawings) will be returned to the author for redraw- 
ing, thus delaying publication of the paper. Photographs 
should be very distinct and show clear black and white con- 
trasts. They must be on glossy white paper. Avoid round and 
oval photographs, 


CAPTIONS AND LEGENDS: Legends or captions must accom- 
pany each drawing or photograph submitted. If written on the 
drawing or photograph, they should be placed below and well out- 
side the part to be reproduced. It is better to place them on 
separate sheets of paper pasted to the back of the drawings or 
photographs. Each table should have a caption such as Table 
1, Table 2, Table 3, etc. Captions should be complete in them- 
selves so as to make the data intelligible to the reader without 
reference to the text. A duplicate list of captious for figures 
should be included as the last page of the manuscript. Use 
“Fig. 1” (not Figure 1), Figs. 3 and 4, etc., in both the text 
and the numbering of illustrations. In the text, ‘‘Eq. (1),”’ or 
“Eqs. (1) and (2)” are preferable to ‘‘Equation (1). In cap- 
tions and legends, except for ‘“‘Fig.” and “‘Eq.”’ and table head- 
ings, write all words in full; do not abbreviate. Avoid placing 
= roe gee written matter in the drawings; it should be in 

e text. 


MATHEMATICAL WorK: Only the very simplest formulae 
should be typewritten; all others should be very carefully written 
in pen and ink, the writing to be large enough so that ample room 
is provided to mark mathematical matter for the printer. A 
considerable space for marking should be allowed above and be- 
low all equations. All complicated equations should be repeated 
on separate sheets with plenty of space left for marking. The 
solidus should be used for simple fractions appearing within 
the text. Make all expressions clear to the typesetter. Greek 
letters used in formulae should be clearly designated by name 
on the margin of the manuscript. All symbols should be clearly 
written and carefully checked. The difference between capital 
and lower-case letters should be clearly distinguished and care 
taken to avoid confusion between zero (0) and the letter (0), 
between the numeral (one) and the letter (ell) and the prime 
(’), between alpha and a, kappa and k, u and mu, v and nu, n 
and eta. All subscripts and exponents should be clearly marked 
and dots and bars over letters or mathematical expressions 
should be avoided. Avoid complicated exponents and subscripts, 
When it is necessary to repeat a complicated expression, it should 
be represented by some convenient symbol. 

NOMENCLATURE AND ABBREVIATIONS: The National Advisory 
Committee for Aeronautics Nomenclature should be used in pref- 
erence to any others. Standard abbreviations should be used, 
and it should be noted that most abbreviations are lower case, 
such as m.p.h., b.m.e.p., i.hp., b.hp., hp., .. . etc. 
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publish and to read of new knowledge and new applications in the field of aeronautical 
engineering. By including in its membership the leading aeronautical specialists in 
other countries an interchange of international thought is made possible. The con- 
tinuous improvement of aircraft depends on scientific research and engineering experi- 
ment applied to design. The Institute desires to bring together all specialists contribut- 
ing to aeronautical progress and, for this purpose, is keeping the dues nominal in amount, 
in view of the fact that practically every specialist already belongs to his own special 
professional society. The dues from members do not provide funds enough both for 
publication of the Journal and for the Institute’s other activities, the difference being 
made up from funds contributed by benefactors and by the aeronautical industry. 


Tue Sxyport, the headquarters of the Institute, located on the fifty-first floor of the 
largest office building in the world, provides members when in New York with a convenient 


and attractive meeting place. 


DUES 
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Members, $10.00; Technical Members, $5.00. 


APPLICATION for MEMBERSHIP 


Application for membership forms may be secured from any member of the Insti- 
tute or from the Secretary. An Applicant is required to furnish the names of several 
members acquainted with his professional work. The Admissions Committee recom- 
mends the grade of membership to the Council which, in its discretion, may extend to the 
Applicant an invitation to join the Institute. 
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